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Abstract 



The International Linear Collider (ILC) is a proposed electron-positron 
collider with a center-of-mass energy of 500 GeV, and a peak luminosity of 
2 ■ 10 34 cm _2 s _1 . The ILC will complement the Large Hadron Collider, a 
proton-proton accelerator, and provide precision measurements, which may 
help in solving some of the fundamental questions at the frontier of scien- 
tific research, such as the origin of mass and the possible existence of new 
principles of nature. 

The linear collider community has set a goal to achieve a precision of 
10~ 4 on the luminosity measurement at the ILC. This may be accomplished 
by constructing a finely granulated calorimeter, which will measure Bhabha 
scattering at small angles. The Bhabha cross-section is theoretically known 
to great precision, yet the rate of Bhabha scattering events, which would 
be measured by the luminosity detector, will be influenced by beam-beam 
effects, and by the inherent energy spread of the collider. The electroweak 
radiative effects can be calculated to high precision and partially checked 
with events with final state photon radiation by distinguishing between the 
observable energy deposits of electrons and of photons in the luminosity 
calorimeter, using a clustering algorithm. 

In order to achieve the design goal, the geometrical parameters of the 
calorimeter need to be reevaluated. This must be done in a generalized 
manner, so as to facilitate future modifications, the need for which is foreseen, 
due to expected changes in the detector concept. 

This work demonstrates that the clustering approach is viable, and that a 
luminosity calorimeter may be designed to match the precision requirements 
on the luminosity measurement. 
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Chapter 1 



Introduction 



1.1 The International Linear Collider 
1.1.1 Physics Case of the ILC 

The triumph of 20 th century particle physics was the development of the 
Standard Model. Experiments determined the particle constituents of ordi- 
nary matter, and identified four forces binding matter and transforming it 
from one form to another. This success leads particle physicists to address 
even more fundamental questions, such as the existence of undiscovered prin- 
ciples of nature, the unification of the four forces, the nature of dark matter 
and dark energy, and the possible existence of extra dimensions. The In- 
ternational Linear Collider (ILC) [1] is expected to play a central role in an 
era of revolutionary advances with breakthrough impact on many of these 
fundamental questions [2]. 

The Standard Model includes a third component beyond particles and 
forces that has not yet been verified, the Higgs mechanism that gives mass to 
the particles. Many scientific opportunities for the ILC involve the Higgs par- 
ticle and related new phenomena at Terascale energies. The Standard Model 
Higgs field permeates the universe, giving mass to elementary particles, and 
breaking a fundamental electroweak force into two, the electromagnetic and 
the weak forces, as is illustrated in Fig. 1.1. Due to quantum effects, the 
Higgs cannot be stable at Terascale energies. As a result, new phenomena, 
such as a new principle of nature called supersymmetry, or extra space di- 
mensions, must be introduced. If the Higgs exists, its properties need to be 
investigated, and if not, then the mechanism responsible for the breaking of 
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the electroweak force must be explained. 
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Figure 1.1: The electromagnetic and weak nuclear forces unify at the Teras- 
cale. 

During the next few years, experiments at the CERN Large Hadron Col- 
lider (LHC), a 14 TeV proton-proton collider, will have the first direct look 
at Terascale physics. If there is a Higgs boson, it is almost certain to be 
found at the LHC and its mass measured by the ATLAS and CMS experi- 
ments [3, 4]. If there is a multiplet of Higgs bosons, there is a good chance 
the LHC experiments will see more than one. However, it will be difficult for 
the LHC to measure the spin and parity of the Higgs particle, and thus to 
establish its essential nature. On the other hand, the ILC, an e + e~ collider, 
can make these measurements accurately, due to the point-like structure of 
the beam particles. If there is more than one decay channel of the Higgs, 
the LHC experiments will determine the ratio of branching fractions (with 
an accuracy roughly of 7 — 30%); the ILC will measure these couplings to 
quarks and vector bosons at the few percent level, and thus reveal whether 
the Higgs is the simple Standard Model object, or something more complex. 

If supersymmetry is responsible for stabilizing the electroweak unification 
at the Terascale and for providing a light Higgs boson, signals of superpartner 
particles should be seen at the LHC. The task would then be to deduce the 
properties of this force, its origins, its relation to the other forces in a unified 
framework, and its role in the earliest moments of the universe. This will 
require precise measurements of the superpartner particles and of the Higgs 



1.1 The International Linear Collider 



3 



particles, and will require the best possible results from the LHC and the 
ILC in a combined analysis. 

Alternative possible structures of the new physics include phenomena con- 
taining extra dimensions, introducing connections between Terascale physics 
and gravity. One possibility is that the weakness of gravity could be un- 
derstood by the escape of gravitons into the new large extra dimensions. 
Events with unbalanced momentum caused by the escaping gravitons could 
be seen at both the LHC and the ILC. Additionally, the ILC could con- 
firm this scenario by observing anomalous electron-positron pair production 
caused by graviton exchange. The advantage of the LHC is a large energy 
coverage. However, in all cases, the ILC is a better tool for understanding 
the background to new physics. 

1.1.2 Design of the ILC 

The ILC is designed to achieve the specifications listed in the ILCSC Pa- 
rameter Subcommittee Report [5]. The three most important requirements 
are 

- an initial center-of-mass energy, y/s, of up to 500 GeV, with the ability 
to upgrade to 1 TeV , 

- a peak luminosity of 2 • 10 34 cm~ 2 s _1 with 75% availability, resulting 
in an integrated luminosity in the first four years of 500 fin" 1 at A /i = 
500 GeV, or equivalent at lower energies, and 

- the ability to scan the energy range 200 < \fs < 500 GeV. 

Additional physics requirements are electron beam polarization > 80%, an 
energy stability and precision < 0.1%, an option for ~ 60% positron beam 
polarization, and alternative e~e~ and 77 collisions. 

The ILC Reference Design Report [6] describes a collider that is intended 
to meet these requirements. Figure 1.2 shows its schematic design. In this 
concept the ILC is based on 1.3 GHz superconducting radio- frequency cav- 
ities operating at a gradient of 31.5 MV/m [7]. The collider operates at a 
repetition rate of 5 Hz with a beam pulse length of roughly 1 msec. The site 
length is 31 km for -y/s = 500 GeV, and would have to be extended to reach 
1 TeV. The beams are prepared in low energy damping rings that operate 
at 5 GeV and are 6.7 km in circumference. They are then accelerated in the 
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main linacs which are ~ 11 km per side. Finally, the beams are focused down 
to very small spot sizes at the collision point with a beam delivery system 
that is ~ 2.2 km per side. 

1.2 The Detector Concept 

Presently, two detectors are considered, with the pull-and-push scheme. With- 
in the next half a year, letters of intent are expected. The European high 
energy committee, which have been working on the fifth version (v5) of the 
so-called "Large Detector Concept" (LDC) [8], has recently joined forces with 
the Japanese and American communities to promote the International Large 
Detector (ILD) [9] concept. Central to the detector design are a micro- vertex 
detector, a time projection chamber (TPC) tracking device, and calorimetry 
with very fine granularity, in order to reconstruct the particle flow for best 
jet-energy resolution. A schematic representation of the LDC(v5) detector 
concept is presented in Fig. 1.3. 

The Forward Region of the ILC 

The instrumentation of the forward region aims at measuring the luminosity, 
providing electron veto at small angles, as well as a fast beam monitoring 
system, and complementing the hermeticity of the full detector. 

The following sub-systems comprise the forward region of the ILC de- 
tector: the luminosity detector (LumiCal) for precise measurement of the 
Bhabha event rate; the beam calorimeter (BeamCal) and the beamstrahlung 
photons monitor (GamCal), for providing a fast feed-back in tuning the lu- 
minosity. BeamCal is also intended to support the determination of beam 
parameters. Both LumiCal and BeamCal extend the angular coverage of 
the electromagnetic calorimeter to small polar angles [10]. The layout of the 
forward region is depicted in Fig. 1.4. 

The requirement for LumiCal is to enable a measurement of the integrated 
luminosity with a relative precision of about 10~ 4 . The use of Bhabha scat- 
tering as the gauge process is motivated by the fact that the cross-section 
is large and dominated by electromagnetic processes, and thus can be calcu- 
lated with very high precision. The purpose of the BeamCal is to efficiently 
detect high energy electrons and photons produced e.g. in low transverse 
momentum QED processes, such as Bhabha scattering and photon-photon 
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Figure 1.2: Schematic design of the ILC for a center of mass energy a/s = 
500 GeV. 
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Figure 1.4: Schematic design of the forward region of the ILC in the LDC(v5) 
detector concept. 
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events. BeamCal is important in order to suppress this dominant background 
in many searches for new particles predicted in scenarios for physics beyond 
the Standard Model. In the polar angle range covered by the BeamCal, typ- 
ically 5 to 45 mrad, high energy electrons must be detected on top of wider 
spread depositions of low energy e + e~ pairs, originating from beamstrahlung 
photon conversions. The measurement of the total energy deposited by these 
pairs, bunch by bunch, can be used to monitor the variation in luminosity 
and provide a fast feedback to the beam delivery system. Moreover, the anal- 
ysis of the shape of the energy flow can be used to extract the parameters 
of the colliding beams. This information can be further used to optimize the 
machine operation. GamCal is used to analyze beamstrahlung photons. It 
will be positioned at a distance of about 180 m from the interaction point. It 
will be sensitive to the energy of the beamstrahlung photon and to the size 
of the beamstrahlung photon cone, which in turn is sensitive to the beam 
parameters. 

1.3 Work Scope 

The focus of this thesis is the design and performance of the luminosity 
calorimeter. The objective is to demonstrate that it is possible to design a 
LumiCal, such that the relative error on the luminosity measurement meets 
the performance requirements. The way to accomplish this is twofold. 

On the one hand, it must be proven that the Bhabha process, which is 
the benchmark process for measuring luminosity at the ILC, may be mea- 
sured directly. This is important for several reasons. For one, the required 
precision with which the Bhabha process needs to be know, is in par with the 
current theoretical uncertainty. Further more, one has to take into account 
that the energy of the colliding beams is not monochromatic. This is mainly 
due to beamstrahlung radiation, energy loss by the incoming positron (elec- 
tron) due to its interaction with the electron (positron) bunch moving in the 
opposite direction. In addition, the physical Bhabha cross-section is affected 
by electroweak radiative effects, which have been calculated to a precision of 
5.3 • 10~ 4 on the Z resonance [11]. A measurement of the differential Bhabha 
cross-section itself would serve as a control mechanism for checking the calcu- 
lations in the region of high momentum photon emission, which is subject to 
the largest corrections. This is done by way of performing clustering in Lu- 
miCal, and measuring the distribution of the (final state) radiative Bhabha 
photons, on top of the electron distribution. 
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The second requirement that needs to be addressed is the measurement of 
the integrated luminosity in LumiCal. The design of LumiCal must balance 
between oftentimes contradicting constraints. On one hand, one would like 
to improve the precision of the luminosity measurement. On the other hand, 
other considerations need to be taken into account, such as minimizing the 
material budget, ensuring the viability of the readout etc. Due to the fact 
that R&D efforts are continuing, the detector concept is still fluid. Conse- 
quently, the restrictions on the size and positioning of LumiCal may change 
every few months. It is, therefore, necessary to define a clear procedure, that 
will allow for adjustment of the design parameters of LumiCal, while keeping 
its performance within the requirements. A study is presented here, in which 
LumiCal is optimized to this effect. The purpose is both to arrive at the 
best design of the calorimeter, under the present constraints imposed by the 
detector concept, and to show how such an optimization should take place. 

This work includes a theoretical introduction (chapter 2 and 3), an overview 
of the clustering algorithm which has been developed for LumiCal, and of its 
performance (chapter 4) 1 , an optimization study of the design of LumiCal 
(chapter 5), and finally, a summary of the results (chapter 6). 



A full description of the clustering algorithm is presented in (Appendix A). 



Chapter 2 
Luminosity 



To measure the cross-section, a, of a certain process we count the number of 
events, N, registered in the detector, and obtain a using the corresponding 
integrated luminosity, L, according to the relation 



Neglecting other systematic uncertainties, the required precision on the lu- 
minosity measurement is given by the statistics of the highest cross-section 
processes which is measured. 

2.1 Luminosity Measurement at the ILC 

2.1.1 Precision Requirements on the Luminosity Mea- 
surement 

At y^i = 340 GeV the cross-section for e + e~ — > W + W~ is about 10 pb, and 
the one for fermion pairs, e + e~ — > q + q~ , is about 5 pb, both scaling with 1/s. 
In both processes one, therefore, expects event samples of O(10 6 ) events in a 
few years of running, which would require a luminosity precision, AL/L, of 
better than 10~ 3 (see Eq. (2.10) below). 

The GigaZ program requires running the collider at an energy correspond- 
ing to the Z pole. The ILC is designed to reach very high luminosity in this 
mode of operation, and thus can become a very powerful laboratory for ad- 
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vancing the tests of the SM, which have been performed at LEP/SLC, to 
a new level of accuracy. The goal of the GigaZ run is a test of the radia- 
tive corrections to the Z-fermion couplings with extremely high precision. 
In general these radiative corrections can be parametrized in terms of three 
parameters, £1,2,3 [12]. Fig. 2.1a shows the expected precision on e lj3 under 
different assumptions [13]. These two parameters can be obtained from the 
Z-observables alone while 62 needs in addition a measurement of the H^-mass. 
Another task at GigaZ is the measurement of the strong coupling constant, 
a s , which can be obtained from the ratio of hadronic to leptonic Z decays 
to a precision of 5 — 7 • 10 -4 . Tests of grand unification are limited by the 
knowledge of the strong coupling constant, as depicted in Fig. 2.1b. Some 
models, e.g. within string theory, predict small deviations from unification, 
thus making this measurement very important. GigaZ is also especially inter- 
esting if no direct evidence for physics beyond the Standard Model is found. 
In this case the structure of radiative corrections should be tested without 
artificial constraints. 




£,*io 3 Q [GeV] 

(a) (b) 

Figure 2.1: (a) Presently known precision for e 13 , expected after the LHC 
and after the ILC measurements. All curves, apart from the one denoted 
"LC, no niw", assume that £2 is equal to its SM value, (b) Unification of 
couplings now and after GigaZ. 

In the GigaZ mode, more than 10 9 hadronic Z decays are expected, which 
would in principle require a luminosity precision of roughly 10~ 5 . However 
there are other systematic uncertainties that play a role. These include the 
selection efficiency for hadronic events and the modification of the cross- 
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section on top of the Breit-Wigner resonance, due to the beam energy spread. 
Hence, a luminosity precision of AL/L ~ 10~ 4 seems adequate [14]. 

2.1.2 Bhabha Scattering as the Gauge Process 

Fig. 2.2 shows the elastic Bhabha scattering process, e + e~ —> e + e~ . 




(a) s channel (b) t channel 



Figure 2.2: Feynman diagrams of the s- and t-channel Born-level elastic 
Bhabha scattering, as denoted in the figure. 

Strictly speaking, Born-level elastic Bhabha scattering never occurs. In 
practice, the process is always accompanied by the emission of electromag- 
netic radiation, for example 

e + e~ e + e"7. (2.2) 

In a simplified picture, a Bhabha event may be depicted as occurring in 
three steps: emission of radiation from the initial particles, Bhabha scat- 
tering, and emission of radiation from the final particles. In the angular 
scattering range considered for the luminosity measurement, one can discard 
the effects of interference between the initial and final state radiation. It 
should also be noted that the initial state radiation is mostly emitted in the 
direction of the beams and travels through the beampipe, thus remaining 
undetected. The ability to distinguish between a final state radiative photon 
and its accompanying lepton is determined by the resolving capabilities of 
the detector, and is a function of the angular separation between the two 
particles. When the two can be resolved, then the experimental measure- 
ments can be compared with the theoretical prediction, and thus the theory 
can be partly tested. 
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The Bhabha scattering includes at the Born level 7 and Z exchange, 
both in the s- and the t-channels. The process may be written in terms of 
ten contributions, where four terms correspond to pure 7 and Z° exchanges 
in the s- and the t-channels, and the other six correspond to 7 — Z° and 
s — t interferences [15]. Taking into account 7 exchange only, the Bhabha 
cross-section may be presented as the sum of three terms, 



dcr B 

In 



a 
2s 



1 + cos 4 (#/2) cos 4 (fl/2) 



1 + cos 2 e 



sin 4 (#/2) 



sin 2 (fl/2) 



(2.3) 



where the scattering angle, 9, is the angle of the scattered lepton with respect 
to the beam, a is the fine structure constant, and s is the center-of-mass en- 
ergy squared. The first and last terms correspond to 7 exchange in the t- and 
s-channels, respectively, and the middle term reflects the s — t interference. 

For small angles (6 < 10°), Bhabha scattering is dominated by the t- 
channel exchange of a photon. Discounting the s-channel contributions, one 
can rewrite Eq. (2.3) in terms of the scattering angle as: 



c?o"b 2ira 2 sin# 327rct 2 1 

~d6 = s sin 4 (#/2) ~ s ¥' ^ 2A > 

The advantage of Bhabha scattering as a luminosity gauge process is 
that the event rate exceeds by far the rates of other physical processes, and 
in addition, the theoretical calculations of the cross-section are under control. 

For the determination of the luminosity, the precise calculation of the 
Bhabha cross-section at small polar angles is needed. Theorists are working 
currently in several laboratories to improve the accuracy of higher-order elec- 
troweak corrections to the Bhabha cross-section [16, 17, 18, 19]. The current 
theoretical uncertainty was estimated to be 5.3- 10~ 4 on the Z resonance [11], 
with the prospect of reducing this uncertainty to 2 ■ 10~ 4 , matching the need 
of GigaZ. 



2.1.3 Beam-Beam Effects at the ILC 

The theoretical uncertainties quoted above are based on the assumption that 
the energy of the colliding beams is monochromatic. For the ILC, though, 
this is not the case. The colliding electron and positron bunches at the ILC 
disrupt one another [20]. Prior to the Bhabha scattering, the interacting 
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particles are likely to have been deflected by the space charge of the opposite 
bunch, and their energies reduced due to the emission of beamstrahlung. To 
take into account the cross-section dependence on s, the probability used 
to produce Bhabha scattering events during the beam-beam collision should 
be rescaled by s/s', where s' is the effective center-of-mass energy after the 
emission of beamstrahlung. The variance in s will, in addition, be aggravated 
by the inherent energy spread of the collider. In general, the collision param- 
eters, such as the size of the collision region and the bunch current, that lead 
to the highest luminosity, also lead to the largest smearing of the luminosity 
spectrum, dL/d\fs. Additionally, the energy measurements can be tempered 
by the presence of beam related backgrounds, such as synchrotron radiation 
and thermal photons of the residual gas, backscattered off the electron beam. 

The acollinearity angle for the e + e~ final state, defined as 



is depicted in Fig. 2.3. Beamstrahlung emissions often occur asymmetrically, 
with either the electron or the positron loosing most of the energy. Hence the 
acollinearity of the final state can be significantly enhanced. The final state 
particles scattered in the acceptance range of LumiCal, following a Bhabha 
interaction, can typically cross a significant part of the opposite bunch. They 
can thus be focused by the electromagnetic field from the corresponding space 
charge, which causes the scattering angle to change. 



e A = e e - + e e+ - n, 



(2.5) 




e 



P. 



e 



Figure 2.3: Definition of the acollinearity angle, #a- 



Both beamstrahlung emissions and electromagnetic deflections vary with 
the bunch length, the horizontal bunch size, and the energy of the collision, 
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and hence so do the resulting biases on the integrated luminosity. Recon- 
structing dL/dy/s from the scattered Bhabha angles is possible [21, 22]. This 
is done by measuring the acollinearity angle, which is related to the difference 
in the energies of the electron and positron beams, in the case of small en- 
ergy and small scattering angle differences. The luminosity spectrum needs 
to be unfolded from the rates for the observed signal- channels in order to pro- 
duce cross-sections as a function of energy. This is especially important for 
such analyses as top-quark and W -boson mass measurements [23]. Knowing 
dL I dy/s also provides a good way to measure the amount of beamstrahlung, 
and thus to predict the corresponding contribution to the bias. 

Contrary to the case with beamstrahlung, there is no direct way to control 
experimentally the bias from the electromagnetic deflections, and so these 
have to be simulated in order to compensate for their effect. 

Since both the beam-beam effects and the collider energy spread depend 
on the parameters of the collisions, it would be very productive to measure 
the Bhabha cross-section itself, and thus better control the systematic errors. 



2.1.4 Relative Error of the Luminosity Measurement 

Several sources (defined below) contribute to the final error of the luminosity 
measurement, 

AL_/AL\ & f AL \ (2 6) 

& V ^ J rec V & J stat 

The typical signature of Bhabha scattering events is the exclusive presence 
of an electron and positron, back to back in the detector. A set of topological 
cuts is applied by comparing the scattering angles of the electron and of the 
positron, and by constraining the difference between, and the magnitude of, 
the energy which is collected in each detector arm [24, 25]. The different 
contributions to the relative error of the luminosity measurement come down 
to 



AL AN 



N - N 

ly rec 1 ' gen 



N 



N, 



yen 



(2.7) 



where N rec and N gen are respectively the number of reconstructed and gen- 
erated Bhabha events, and 9 m i n and 9 max are the respective low and high 
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bounds on the fiducial volume (acceptance region) of the detector. 

An error on the measurement is, therefore, incurred when events are 
miscounted. This may happen for several reasons. One of the causes for 
miscounting events has to do with knowledge of the effectiveness of the cuts 
for distinguishing between Bhabha events and background; the efficiency 
and purity of the cuts must be known to good precision in order to avoid 
counting errors. Another important factor is migration of events out of the 
acceptance region and into it. This may occur if a shower's position in 
the detector is not reconstructed well, resulting in inclusion of events which 
were in actuality out of the fiducial volume, or visa-versa. Errors can also 
result from poor knowledge of the geometrical properties of the detector. For 
instance, displacement of the two arms of the detector with respect to each 
other, or with regard to the interaction point, may lead to systematic biases 
in the position reconstruction. 

In the following, the errors which were discussed above are quantified. 



Error in Reconstruction of the Polar Scattering Angle 

The polar angle dependence of the Bhabha cross-section is da/d0 ~ 9~ 3 
(Eq. (2.4)). This means that the total Bhabha cross-section within the an- 
gular range [9 min ,9 max ] is 



where the 6 max dependence can be neglected. The relative error on luminosity 
is proportional to the relative error on the Bhabha cross-section, 



fAL\ A9 , . 

(it) =2 ^- (2 - 9) 

The analytic approximation of Eq. (2.9) has been shown to hold well in 
practice [25]. Its implication is that the polar bias, A9, and the minimal 
polar bound of the fiducial volume, 9 m i n , are the two most important param- 
eters that affect the precision of the luminosity measurement. The steep fall 
of the Bhabha cross-section with the polar angle translates into significant 
differences in the counting rates of Bhabha events, for small changes in the 
angular acceptance range. 
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Statistical Error of the Number of Expected Bhabha Events 

The probability of observing Bhabha scattering in a given event is determined 
by the Poisson distribution. The variance of the distribution is then equal 
to the average number of observed Bhabha scatterings, N. The relative 
statistical error stemming from Eq. (2.1), is, therefore 



Equation (2.10) is the driving force behind the precision requirements, 
which were stated in Sect. 2.1.1. 

Additional Sources of Error 

In congruence to the two major sources of error which were discussed above, 
several other factors need to be considered in order to keep the design goal of 
AL/L ~ 10 -4 . These include controlling the position of the inner radius of 
the detector on a /im precision level, and the distance between its two arms 
to within 0(100 /im), to name just two. A full account may be found in [24]. 

2.2 The Luminosity Calorimeter 

LumiCal is a tungsten-silicon sandwich calorimeter. In the present ILD lay- 
out the detector is placed 2.27 m from the interaction point. The LumiCal 
inner radius is 80 mm, and its outer radius is 350 mm, so that its polar 
coverage is 35 to 153 mrad. The longitudinal part of the detector consists of 
layers, each composed of 3.5 mm of tungsten, which is equivalent to 1 radia- 
tion length thickness. Behind each tungsten layer there is a 0.6 mm ceramic 
support, a 0.3 mm silicon sensors plane, and a 0.1 mm gap for electron- 
ics. LumiCal is comprised of 30 longitudinal layers. The transverse plane 
is subdivided in the radial and azimuthal directions. The number of radial 
divisions is 104, and the number of azimuthal divisions is 96. Figure 2.4 
presents the segmentation scheme of LumiCal. 

The two half barrels can be clamped on the closed beam pipe. The posi- 
tion of the two parts of the detector with respect to each other will be fixed by 
the help of precise pins placed at the top and bottom of each C shaped steel 
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Figure 2.4: Half plane of LumiCal silicon sensors (every fourth raidal segment 
is drawn). 



frame [26]. The latter stabilizes the structure and carries the heavy tung- 
sten disks by the bolts. The gravitation sag of the tungsten absorber can be 
kept in required tolerance [27]. The silicon sensors are glued to the tungsten 
surface with capton foil insulation. Space for readout electronics, connectors 
and cooling is foreseen at the outer radius of the calorimeter. The sensor 
plane will be built from a few tiles because the current technology is based 
on 6-inch wafers, and at the moment it is unclear if and when larger wafers 
will be available. The tiles of the silicon sensors will be glued to a thick film 
support ceramic plate or directly to a tungsten surface with some insulation. 
Reference marks are foreseen on the detector surface for precision position- 
ing. The layout of the sensors and the mechanical design of the calorimeter 
does not allow for sensors to overlap. To reduce the impact of the gaps, odd 
and even planes are rotated by half the azimuthal cell pitch (1.875°). The 
silicon diodes will be usual planar high resistivity silicon sensors. Figure 2.5 
presents the foreseen mechanical design of LumiCal. 
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Figure 2.5: Foreseen mechanical design of LumiCal. 



Chapter 3 

Development of EM Showers in 
LumiCal 

3.1 Basic Concepts in Calorimetry 

3.1.1 Energy Loss by Electrons and Photons 

High-energy electrons predominantly loose energy in matter by bremsstrahlung, 
and high-energy photons by e + e~ pair production. The characteristic amount 
of matter traversed for these related interactions is called the radiation length, 
X . It is both the mean distance over which a high-energy electron looses all 
but 1/e of its energy by bremsstrahlung, and | of the mean free path for pair 
production by a high-energy photon [28]. The radiation length is also the 
appropriate scale length for describing high-energy electromagnetic showers. 

At low energies, electrons and positrons primarily loose energy by ion- 
ization, although other processes (M0ller scattering, Bhabha scattering, e + 
annihilation) contribute as well, as shown in Fig. 3.1a. While ionization 
loss-rates rise logarithmically with energy, bremsstrahlung losses rise nearly 
linearly (fractional loss is nearly independent of energy), and dominate above 
a few tens of MeV in most materials. 

At low energies, the photon cross-section is dominated by the photoelec- 
tric effect, although Compton scattering, Rayleigh scattering, and photonu- 
clear absorption also contribute. The photoelectric cross-section is character- 
ized by discontinuities (absorption edges) as thresholds for photoionization 
of various atomic levels are reached. The increasing dominance of pair pro- 
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duction as the energy increases is shown in Fig. 3.1b. The cross-section is 
very closely related to that for bremsstrahlung, since the Feynman diagrams 
are variants of one another. 

3.1.2 Electromagnetic Showers 

When a high-energy electron or photon is incident on a thick absorber, it initi- 
ates an electromagnetic (EM) shower as pair production and bremsstrahlung 
generate more electrons and photons with lower energy. The longitudinal de- 
velopment is governed by the high-energy part of the cascade, and therefore 
scales as the radiation length in the material. Electron energies eventually 
fall below the critical energy (defined below), and then dissipate their energy 
by ionization and excitation, rather than by the generation of more shower 
particles. 

The transverse development of electromagnetic showers scales fairly ac- 
curately with the Moliere radius, Rm, given by [29] 

Rja = X — , (3.1) 

where E s ~ 21 MeV, and E c is the critical energy, which is defined as the 
energy at which the ionization loss per radiation length is equal to the electron 
energy [30]. On average, only 10% of the energy of an EM shower lies outside 
a cylinder with radius Rm around the shower-center. 

3.2 Simulation of the Detector Response 

The response of LumiCal to the passage of particles was simulated using 
Mokka, version 06-05-p02 [31]. Mokka is an application of a general pur- 
pose detector simulation package, GEANT4, of which version 9.0.p01 was 
used [32]. The Mokka model chosen was LDC00_03Rp, where LumiCal is 
constructed by the LumiCalX super driver. The output of Mokka is in the 
LCIO format. Several Marlin processors were written in order to analyze 
the LCIO output. Marlin is a C++ software framework for the ILC soft- 
ware [33]. It uses the LCIO data model and can be used for all tasks that 
involve processing of LCIO files, e.g. reconstruction and analysis. The idea 
is that every computing task is implemented as a processor (module) that 
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Figure 3.1: (a) Fractional energy loss per radiation length in lead as a func- 
tion of electron or positron energy. Electron (positron) scattering is consid- 
ered as ionization when the energy loss per collision is below 0.255 MeV, and 
as M0ller (Bhabha) scattering when it is above, (b) Probability P that a 
photon interaction will result in conversion to an e + e~ pair (the figures are 
taken from [29]). 
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analyses data in an LCEvent and creates additional output collections that 
are added to the event. Version 00-09-08 of the program was used. 

The geometry of LumiCal which was simulated is that which is described 
in Sect. 2.2. 

Figure 3.2 shows the generation spectrum of electrons, positrons and pho- 
tons for a 250 GeV EM shower in LumiCal 1 . These particles traverse the 
layers of tungsten and deposit energy in the silicon sensors mainly through 
ionization. In Fig. 3.2a each entry represents the z-position (relative to the 
IP) of the creation of a shower particle with a given energy. In Fig. 3.2b a 
normalized profile of the energy as a function of the distance is presented. 
As the shower develops in depth in the calorimeter, new shower particles are 
created with less and less energy. Eventually the energy falls off below the 
threshold of ionization. 

Two normalized distributions are overlaid in Fig. 3.3a, the number of 
shower particles and the deposited energy, both as a function of the layer 
number in LumiCal. Electron showers of 250 GeV were used. The energy 
deposited in the silicon sensors is proportional to the number of charged 
shower particles. This is consistent with the fact that both distributions 
have a similar shape. However, while the distribution of the number of 
shower particles peaks at the ninth layer, that of the energy deposition peaks 
at the tenth. A displacement between the distributions by one layer, which is 
equivalent to one radiation length, is apparent. This is due to the fact that 
part of the shower is comprised of photons, which do not deposit energy, 
but are later converted to electron-positron pairs, which do. In Fig. 3.3b is 
presented the normalized distribution of the number of cell hits for 250 GeV 
electron showers as a function of the layer in LumiCal. The number of cells 
which register a hit peaks around layer number 13. 

Figure 3.4 shows the profile of the energy deposited in LumiCal for a 
single 250 GeV electron shower. Integration is made along the ^-direction 
for equivalent values of the x and y coordinates, which are taken as the 
centers of the cells of the relevant hits. The global shower- center is defined 
as the center of gravity of the shower profile, using cell energies as weights. 
The polar symmetry of the segmentation of the detector is also visible in the 
figure, and it is evident that LumiCal is more finely granulated in the radial 
direction than in the azimuthal direction. 

1 The shower also contains protons and neutrons. The contribution of these particles 
to the energy deposited in LumiCal is negligible due to their relative low number, and so 
will not be discussed here. 
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(a) (b) 

Figure 3.2: (a) Distribution of the creation position relative to the IP and of 
the energy of shower particles in LumiCal for 250 GeV electron showers, (b) 
Normalized distribution of the energy at the creation vertex of shower parti- 
cles, as a function of the point of creation. The creation point is represented 
by corresponding scales of the distance from the IP and the LumiCal layer 
number. 
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Figure 3.3: (a) Normalized distributions of the number of shower particles 
and of the energy deposited in the silicon sensors of LumiCal as a function 
of the layer, as denoted in the figure. Electron showers of 250 GeV were 
simulated, (b) Normalized distribution of the number of cell-hits for 250 GeV 
electrons showers as a function of the layer in LumiCal. 
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Figure 3.4: Transverse energy profile for a 250 GeV electron shower in Lu- 
miCal. The red circle represents an area bound within one Moliere radius 
around the shower center. 



Figure 3.5a shows the distribution of the distance around the global 
shower-center, in which 90% of the integrated shower energy may be found. 
The distribution is centered around 14 mm, which is, by definition, the 
Moliere radius of LumiCal, R^. Taking into account all of the hits of the 
shower, the polar and azimuthal production angles of the initiating particle 
may be reconstructed. Local shower-centers are defined on a layer-to-layer 
basis as the extrapolation of the trajectory of the particle according to these 
angles. We define the distance around the local shower-center of layer £, 
in which 90% of the layer's energy is deposited, as its layer-radius, r(£). 
Figure 3.5b shows the dependence of r(£) on the layer number, I. 

According to the distributions in Fig. 3.3, in the first layers there are 
few cells which register hits. For this reason there is no clear local shower- 
center, and the area that encompasses 90% of the energy of the layer is large. 
The information in these layers is, therefore, not sufficient to obtain a clear 
description of the shower. This effect is lessened as the shower develops in 
depth and the number of cell-hits increases. Starting at the fifth layer, the 
shower becomes homogeneous. Beyond this point the shower becomes more 
and more wide spread in depth, and its diameter may be estimated to good 
approximation by a power-law. For layer numbers higher than 16, the shower 
exceeds 150% of Rm and looses homogeneity once again. This behavior is 
supported by Fig. 3.3, which shows that for layers beyond the shower-peak, 
the number of shower-particles falls off faster than the number of cells which 
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Figure 3.5: (a) Distribution of Rm, the distance around the global-shower 
center, in which 90% of the integrated shower energy may be found, (b) 
Dependence on the layer number, £, of the layer-radius, r(£), which is the 
distance around the local shower-center in which 90% of the energy of a layer 
is deposited. 
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are hit. Since the shower becomes attenuated for high layer numbers, it 
is difficult to determine the local shower-center with good accuracy. It is, 
therefore, useful to define an effective layer-radius, r e ff(£), by extrapolation 
of the behavior of r(l) at the middle layers, to the front and back layers, 

r eff {£) = r(6 < I < 24) « 5 + 0.23 • I + 0.04 ■ I 2 (3.2) 

In summary, the transverse profile of EM showers in LumiCal is char- 
acterized by a peak of the shower around the tenth layer. The number of 
shower-particles before layer six and after layer 16 is small compared to that 
in the inner layers. The energy of shower particles degrades in depth, and 
accordingly so do the energy deposits, while the shower becomes more wide 
spread. The front layers (layers one to five) are, therefore, characterized 
by a small number of concentrated energy deposits. The middle, so-called, 
shower-peak layers (layers six to 16) "register large energy contributions, and 
the back layers (layers 16 to 30) are characterized by a decreasing number of 
low-energy shower particles, which deposit little energy in a dispersed man- 
ner. The shower has a prominent center, within Rj^ (14 mm) of which most 
of the shower energy is concentrated. On a layer-by-layer basis, most of the 
energy may be found within an effective layer-radius from the center, which 
is parameterized by Eq. (3.2). 



On an event- by- event basis the longitudinal profile is not always as smooth as the one 
represented in Fig. 3.3a. As a result, the shower-peak layers are not necessarily consecutive. 



Chapter 4 



A Clustering Algorithm for 
LumiCal 

In the running conditions of the accelerator, the selection of Bhabha can- 
didate events will require pattern recognition in the main detector and in 
LumiCal. Here the first attempt to perform clustering in LumiCal is pre- 
sented. The main focus is on clustering optimized for EM showers, with the 
intent of resolving events in which hard photons were emitted in the final 
state. 

As explained in chapter 3, high energy electrons and photons which tra- 
verse LumiCal loose energy in the tungsten layers mainly by the creation of 
electron-positron pairs and of photons, which in turn also loose energy by 
the same processes. The cascade of particles is propagated until most of the 
energy of the initial particle has been absorbed in the calorimeter. These 
secondary particles make up an electromagnetic shower that is sampled in 
the silicon sensors that make up the back side of each layer. When two (or 
more) high energy particles enter LumiCal, an EM shower will develop for 
each particle, and the multiple showers will overlap to some degree, depend- 
ing on the initial creation conditions of each shower. The ability to separate 
any pair of showers is subject to the amount of intermixing of the pair. 

4.1 Outline of the Clustering Algorithm 

The clustering algorithm which was designed for LumiCal was written as a 
series of Marlin (version 00-09-08) processors [33]. It operates in three main 
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phases, 

- selection of shower-peak layers, and two-dimensional clustering therein, 

- fixing of the number of global (three-dimensional) clusters, and collec- 
tion of all hits onto these, 

- characterization of the global-clusters, by means of the evaluation of 
their energy density. 

A short description of each phase will now follow. A complete account 
can be found in Appendix A. 

4.1.1 Clustering in the Shower-Peak Layers 

In the first layers of LumiCal, only a few hits from the shower are registered, 
as was discussed above. In addition to the hits from the main showers, there 
may also be contributions owing to backscattered particles or background 
processes. These particles have low energy and do not propagate to the 
inner layers, but their energy is of the order of the depositions of the showers 
of interest. In order to make a good estimate of the number of main showers, 
one must, therefore, begin by considering the information in the shower-peak 
layers. This process is done in two steps, which are described below, near- 
neighbor clustering and cluster-merging. 

Near-Neighbor Clustering 

Initially, clusters are created from groups of closely-connected cells. This 
is done by means of the method of near-neighbor clustering (NNC), which 
exploits the gradient of energy around local shower-centers. The assumption 
is that in first order, the further a hit is relative to the shower center, the 
lower its energy. By comparing the energy distribution around the center 
at growing distances, one may check whether the energy is increasing or 
decreasing. An increase in energy for growing distance from the shower- 
center would then imply that the hit should be associated with a different 
shower. 

For each shower-peak layer separately, the algorithm associates each cell 
which has an energy deposit with its highest-energy near-neighbor. The 
result of the NNC phase is a collection of clusters in each layer, centered 
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around local maxima, as illustrated for a single layer in Fig. 4.1a. In this 
example the algorithm produces six clusters, which are enumerated in the 
figure. The different clusters are also distinguished by different color groups, 
where darker shadings indicate a higher energy content of the cell in question. 

Cluster-Merging 

The next step in the algorithm is cluster-merging. The NNC method only 
connects cells which are relatively close, while showers tend to spread out over 
a large range of cells, as indicated by Fig. 3.5b. The cluster-merging proce- 
dure begins by assigning a center-position to each existing cluster. Weights 
are then computed for each cluster to merge with the rest of the clusters. 
In general, the weights are proportional to the energy of the candidate clus- 
ter, and inversely proportional to the distance between the pair of clusters. 
Several variations of the weighting process are tested in consecutive merg- 
ing attempts. The result of the algorithm after the cluster-merging phase is 
illustrated for a single layer in Fig. 4.1b. 

4.1.2 Global Clustering 

The most important stage of the clustering algorithm is the determination of 
the number of reconstructed showers. The aftermath of the clustering in the 
shower-peak layers is several collections of two-dimensional hit aggregates, 
the number of which varies from layer to layer. The final number of showers 
is then determined as the most frequent value of the layer-cluster number, 
derived from the collections in the shower-peak layers. 

Once the number of global-showers is fixed, cells from non-shower-peak 
layers are associated with one of the global-showers. This is done by ex- 
trapolating the propagation of each shower through LumiCal, using the in- 
formation from the shower-peak layers. Following the extrapolation, cells 
are merged with the extrapolated global-cluster centers in each layer. This 
process is facilitated by assuming a typical shower-size, defined according 
to the parameterization of Eq. (3.2) (Sect. 3.2), which acts as a temporary 
center-of-gravity. Once the core cells within the assumed shower-radius are 
associated with the global-centers in these layers, the rest of the cells may also 
be added. A weighing method, similar to the one used in the shower-peak 
layers, is used here as well. 

Figure 4.2 shows a schematic representation of the global-clustering phase. 
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(b) 

Figure 4.1: (a) Schematic representation of the results of the NNC phase 
for a single layer. Six clusters are found by the algorithm. The different 
clusters are enumerated, as well as distinguished by different color groups. 
Darker shadings indicate a higher energy content of the cell in question, 
(b) Evolution of the results from (a) after the cluster-merging phase. Two 
clusters remain after clusters two, three and five were merged with cluster 
one, and cluster four was merged with cluster six. 
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LumiCal layers are represented by the large blue disks, and layer-clusters are 
represented by the small triangles, squares and circles. Three layers have two 
layer-clusters, one layer has three layer-clusters and one layers has one layer- 
cluster. The first and last layers have no layer-clusters, since they are not 
shower-peak layers. The global number of clusters is two, and the layer clus- 
ters are associated with each other according to the straight lines. The lines 
also define the global-cluster positions in the non- shower-peak layers. The 
cluster represented by a circle in the layer, in which three layer-clusters were 
found, will be disbanded. Its hits will be associated with either the "square" 
or "triangle" global-clusters. The layer-cluster in the layer, in which only 
one cluster was found, will also be disbanded. The hits will then be clustered 
around the virtual-cluster positions, represented by the intersection of the 
straight lines with the layer. A similar procedure will also be performed in 
the first and last layers, where no layer-clusters were constructed previously. 




Figure 4.2: Schematic representation of the global-clustering phase. LumiCal 
layers are represented by the large blue disks, and layer-clusters are repre- 
sented by the small triangles, squares and circles. The straight lines show 
the extrapolated position of the global-clusters in all layers. 
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4.1.3 Corrections Based on the Energy Distribution 

At this point all of the hits in the calorimeter have been integrated into one 
of the global-clusters. Before moving on, it is beneficial to make sure that 
the clusters have the expected characteristics. 

Energy Density Test 

The EM shower development in LumiCal has been described in Sect. 3.2. 
Accordingly, one would expect that 90% of a cluster's energy would be found 
within one Moliere radius, Rm, of its center. While statistically this is true, 
on a case-by-case basis fluctuations may occur, and thus it should not be set 
as a hard rule. It is possible to define, instead, a set of tests based on the 
amount of energy which is located in proximity to each cluster-center. Both 
the energy-density of individual clusters and that of all the clusters together 
is evaluated. When the global-clusters fall short, a quick re-clustering is 
possible. The first step is to profile the energy of all cells in the longitudinal 
direction (see Fig. 3.4), and then strip away low energy cell contributions and 
perform the profiling procedure again in successive iterations. This way, it is 
sometimes possible to reliably locate the high density shower centers. Global 
clusters are then constructed around these centers. The energy density of 
the new clusters is compared to that of the original clusters, and the best set 
is finally kept. 

Unfolding of Mixed Clusters 

Another modification that can be made in the aftermath of the clustering 
procedure, is allocation of hits for mixed cluster pairs. When a pair of showers 
develop in close proximity to each other (in terms of their Moliere radius), 
some cells receive energy depositions from both showers. The problem is, that 
the clustering procedure associates each cell with only one cluster. This biases 
the energy content, especially of low-energy clusters, due to the fact that their 
energy tends to be greatly over-estimated by contributions from high-energy 
clusters. High-energy clusters are less affected, because percentage-wise, the 
variance in energy caused by low-energy clusters is insignificant. A way to 
correct for this effect is to evaluate the energy distribution of each cluster in 
the region furthest away from the position of its counterpart. If one assumes 
that the shape of each shower is smooth 1 , the distribution of hits in the 
area where there is no mixing can be used to predict the distribution in the 
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mixed area. Correction factors are then derived on a cell-by-cell basis, and 
the energy is split between the pair of clusters accordingly. 

4.2 Physics Sample 

The physics sample which was investigated consisted of 3-10 4 Bhabha scatter- 
ing events with center-of-mass energy y/s = 500 GeV. The events were gen- 
erated using BHWIDE, version 1.04 [34]. BHWIDE is a wide angle Bhabha 
MC, which contains the electro-weak contributions, which are important for 
the high energy e + e~ interactions considered here. The sample contains only 
events in which the leptons are scattered within the polar angular range 
35 < 9 < 153 mrad. While all of these events were processed by the clus- 
tering algorithm, some were eventually discarded. Only events in which the 
reconstructed cluster with the highest energy content was found within the 
fiducial volume (acceptance range) of LumiCal, 41.5 < 9 < 131 mrad, were 
kept. Individual clusters were constrained in the same way. The reason for 
this is that showers whose position is reconstructed outside the fiducial vol- 
ume are not fully contained, i.e. a significant amount of their energy leaks 
out of the detector, making the reconstruction process unreliable. 

Figure 4.3 shows the energy spectrum of the scattered leptons and radia- 
tive photons. The lepton distribution peaks at 250 GeV, as expected, and 
has a long tail of lower energies, accounting for the energy which was carried 
away by the photons. 

Figure 4.4 shows the polar 11 and azimuthal production angles, 9 and <fr, of 
scattered leptons and radiative photons. The distribution of the polar angle 
is cut according to the fiducial volume of LumiCal. As expected in light of 
Eq. (2.4), the distribution of the polar angle falls off rapidly with 9, and the 
distribution of the azimuthal angle is flat. 

Since most initial state radiative photons travel through the beampipe 
and are undetected (see Sect. 2.1.2), only final state photons are considered. 
Conservation of momentum dictates that the more energy these photons take 

I In fact, the assumption of smoothness is not always correct. This is due to statistical 
fluctuations in the shower development, and also to the fact that the difference of the 
cell sizes in play are not taken into account. Despite this, the method does improve the 
estimation of cluster energy 

11 Naturally the electron and the positron have polar angles of opposite signs, but as 
the distributions of the production angles are equivalent for either one, this sign will be 
ignored throughout the following. 
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Figure 4.3: Distribution of the production energy of scattered leptons and 
radiative photons for Bhabha scattering events with center-of-mass energy 
y/s = 500 GeV in the LumiCal fiducial volume. 
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Figure 4.4: Distributions of the polar (a) and azimuthal (b) production an- 
gles, 9 and 0, of leptons and photons, as denoted in the figures. The Bhabha 
scattering events were simulated with center-of-mass energy y/s = 500 GeV 
in the LumiCal fiducial volume. 
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from the lepton, the smaller the angular separation between the two. This 
is confirmed by Fig. 4.5a, which shows the correlation between the photon 
energy and its angular separation from the accompanying lepton, Af^ j7 . In 
Fig. 4.5b the energy dependence of Fig. 4.5a is integrated and normalized, 
showing the event rate for the photon emission. The distance in this case 
is expressed as the separation between the pair of particles on the face of 
LumiCal in units of mm and of Moliere radius. It is apparent from the 
distributions that the vast majority of radiative photons is of low energy, 
and enter LumiCal in close proximity to the lepton. 
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Figure 4.5: (a) Correlation between the angular separation between lep- 
tons and radiative photons, Af^ i7 , and the photon energy. The spectrum 
of particles is generated for Bhabha scattering with center-of-mass energy 
y/s = 500 GeV in the LumiCal fiducial volume, (b) Energy profile of the 
distribution in (a). The distance in this case is expressed in units of mm, as 
the separation between the pair of particles on the face of LumiCal, d er/ . An 
equivalent scale is also shown in units of Moliere radius, Rjy[. 



4.3 Performance of the Clustering Algorithm 

The distributions for the position and energy presented in the previous sec- 
tion were drawn from the raw output of the BHWIDE event generator. As 
such, they represent an ideal description of Bhabha scattering. In reality, 
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observables are distorted by the inherent resolution of the measuring device. 
The energy resolution, which is determined by the amount of leakage, and by 
the sampling rate of the calorimeter (see chapter 5), incurs an error on the 
signal-to-energy calibration of LumiCal. Similarly, the polar and azimuthal 
reconstructed angles have a resolution, and also a bias, of their own. In 
order to analyze the output of the clustering algorithm it is necessary to iso- 
late the errors in reconstruction resulting from the clustering, from the other 
systematic uncertainties of LumiCal. 

To this effect, two classes of objects may be defined. The basic simulation- 
truth data will be represented by showers, which contain all of the hits which 
belong to an EM shower initiated by a single particle. These will be referred 
to as generated showers. Since a single detector cell may contain contribu- 
tions from more than one EM shower, generated showers may share cells. 
Hit collections, built by the clustering algorithm, will be referred to as re- 
constructed clusters. In order to remove the systematic uncertainties, the 
properties of both the showers and the clusters are reconstructed in the same 
manner, using information from the detector cells. 

Since there is no way to distinguish in practice between EM showers ini- 
tiated by leptons and those started by photons, reconstructed clusters and 
generated showers will be referred to as having either high- energy, or low- 
energy, which correspond to effective leptons, and effective photons respec- 
tively. High-energy clusters (showers) are identified as those that have the 
highest integrated energy content among the set of all reconstructed clus- 
ters (generated showers). The rest of the clusters (showers) are identified as 
low-energy clusters (showers). 

4.3.1 Event Selection 

Figure 4.6 shows the success and failure of the clustering algorithm in distin- 
guishing between a pair of generated showers as a function of the separation 
distance between the pair, d pair , and of the energy of the low-energy shower, 
Ei ow . Failure of the algorithm may take two forms. A pair of generated 
showers may be merged into one reconstructed cluster (Fig. 4.6a), or one 
shower may be split into two clusters (Fig. 4.6b). As expected, since the 
great majority of radiative photons enter LumiCal within a small distance 
from the leptons, separation between the showers of the two particles is not 
trivial. The difficulty is enhanced due to the increasing size of showers, as 
they develop in depth in LumiCal. Distinguishing between pairs of showers 
becomes easier when either Ei ow or d pa i T increase in value. 
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Figure 4.6: Success and failure of the clustering algorithm in distinguishing 
between a pair of generated showers as a function of the separation distance 
between the pair, d pa i r , and of the energy of the low-energy shower, E\ ow . 
The distance d pa i r is expressed in units of mm and of Moliere radius, Rm- 
Both figures show success (1 — > 1) of the algorithm; (a) also shows cases 
where a pair of generated showers are merged into one reconstructed cluster 
(2 — > 1), and (b) also shows cases where a single shower is split into two 
clusters (1 — > 2). 
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It is, therefore, required to set low bounds on the energy of a cluster, and 
on the separation between any pair of clusters. When the algorithm produces 
results that do not pass the cuts, the two clusters are integrated into one. 
In order to compare with theory the distribution of clusters after making 
this merging- cut on Ei ow and d pair , one must also apply the same restric- 
tions on the generated showers. The generated showers follow a distribution 
complying with an effective Bhabha cross-section. 

The distinction between the original and the effective cross-sections is 
important, and it must be noted that the effective cross-section can only be 
computed by simulating the detector response. The position of a cluster is 
reconstructed by making a cut on cell energy, relative to the entire cluster en- 
ergy (see Eq. (A. 2)). As a result, an integration of a pair of clusters into one, 
sets the position of the merged cluster to an a-priori unpredictable value 111 . 
The momentum of the initiating particles will, in some cases, not balance 
with that of the effective (merged) particle. Summing up deposits from 
multiple showers in LumiCal is, therefore, not equivalent to any summation 
procedure that might be done on the cross-section, at the generated-particle 
level. 



4.3.2 Observables 
Quantification of the Performance 

The error on the effective cross-section will depend on the number of mis- 
counted showers. In order to judge the success of the algorithm, one may 
evaluate its acceptance, A, purity, CP, and efficiency, £, which are defined as 



■A = Nl ^ AT , CP = AT Nl ^ AT , and £ = (4.1) 

where is the number of generated showers which were reconstructed as 
one cluster by the algorithm, N 2 ^i is the number of pairs of showers which 
were reconstructed as one cluster, and iVi-^ is the number of single showers 
which were separated into two reconstructed clusters. 



111 For instance, if a cluster is of much higher energy than its counterpart, then the energy 
contributions of the low-energy cluster will not be taken into account in the position 
reconstruction at all. For the reconstruction of single showers this is not a problem, since 
the shower is of homogeneous shape around a defined center. For the case of two showers, 
which are far apart, this is no longer the case. 
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The values of the acceptance, purity and efficiency are presented in Ta- 
ble 4.1 for several pairs of merging-cuts on the minimal energy and the sep- 
aration distance between a pair of clusters. Also shown is the fraction of 
radiative photons which are available for reconstruction after applying the 
merging-cuts, 



iV 7 (cut) 
JV 7 (all) 



(4.2) 



where iV 7 (all) is the total number of radiative photons in the fiducial volume 
of LumiCal, and iV 7 (cut) is the number of photons in LumiCal which also 
pass the merging-cuts on Ei ow and d pa i r - 



Cuts 


p 7 [%} 


A [%) 


9 [%] 


£ [%] 


dpair [Rm] 


E min [GeV] 


0.5 


25 


6.6 


69 


96 


71 


0.75 


20 


5.9 


85 


95 


90 


0.75 


25 


5.2 


58 


96 


89 


1 


15 


6 


94 


93 


100 


1 


20 


5.2 


95 


95 


100 


1 


25 


4.6 


95 


96 


98 


1.5 


20 


4.3 


99 


98 


100 



Table 4.1: The values of the percentage of photon showers, which are available 
for reconstruction, p 7 , and of the acceptance, A, purity, CP, and efficiency, £, 
of the algorithm, for several pairs of merging-cuts on the minimal energy of 
a cluster, E m i n , and on the separation distance between a pair of clusters, 
dpair- The merging-cut d pa i r is expressed in units of the Moliere radius, Rm- 

The relative error of the effective cross-section as a result of miscounting 
depends on the observed number of effective leptons and photons, and on 
the fractions of miscounted events out of the relevant event population. The 
probability of finding a given value for N±^2 or for iV^i is given by the 
binomial distribution, and so the relative error on either one is 
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AiV,, 



\/Ni n pq 




(4.3) 



where p is the probability to miscount in a given event, q = 1—p, and Ni a is 
either the number of effective leptons, Ni, or the number of effective photons, 
Nj, depending on the type of miscounting IV . Values for p and q were derived 
from running the clustering algorithm on the sample of Bhabha events with 
different sets of merging-cuts on Ei ow and d pa i r . The corresponding relative 
errors are shown in Table 4.2. Also shown there is the relative error 



which corresponds to the total error resulting from both types of miscounting, 
rescaled for an integrated luminosity of 500 fb _1 . 

It is apparent from Tables 4.1 and 4.2 that achieving a minimum of error 
in counting the number of effective photons, depends both on the size of 
the sample of available photons, and on the sensitivity of the algorithm to 
miscounting. For merging-cuts in energy > 20 GeV and distance > Rm, 
the algorithm makes relatively few mistakes. The decision on where exactly 
to set the merging-cuts reduces to the choice of maximizing the measurable 
amount of statistics. 

Event-by-Event Comparison of Observables 

Other than counting the number of low and high-energy clusters and compar- 
ing the results to the expected numbers, deduced from the effective Bhabha 
cross-section, the properties of the clusters may also be evaluated. For this 
purpose, one may produce such distributions as the production angles of 
clusters, the angular separation between pairs of clusters, and the value of 
cluster-energy. A first step in this process is to look at the shower /cluster 
differences on an event-by-event basis. 

The energy of the particle which initiated a generated shower (recon- 
structed cluster) is determined by integrating all the contributions of the 

IV When computing the number of single showers which were split into two clusters, 
— > Ni, since the candidate for false splitting will come from the population of single 
(lepton) showers. On the other hand, merging of two showers into one may only happen 
when a photon shower exists, so that in this case Ng_ y N 7 . This distinction is important, 
as the number of effective leptons far outweighs the number of effective photons. 




(4.4) 
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Cuts 


AiVi_> a 




AN tot 


dpair [Rm] 


E ml n [GeV] 








0.5 


25 


4.2 • 10" 4 


31.5 • 10- 2 


10.3 ■ 10~ 5 


0.75 


20 


7.6 ■ 10" 4 


14.6 • 10~ 2 


7.5- 10~ 5 


0.75 


25 


5.4- 10~ 4 


14.6 • lO" 2 


8- 10~ 5 


1 


15 


12.9 • 10~ 4 


6.3 ■ 10~ 2 


5.1 • 10~ 5 


1 


20 


7- 10" 4 


4.6 • 10~ 2 


4.9 • 10~ 5 


1 


25 


5.1 ■ 10" 4 


5.2 • lO" 2 


5.3- 10~ 5 


1.5 


20 


3.2 ■ 10" 4 


0.5 ■ 10- 2 


1.8 • 10~ 5 



Table 4.2: The relative errors on the miscounting of clusters (using N^ 2 
and N2-+1), and the total relative error on the measurement of the effective 
Bhabha cross-section (using N to t), for several pairs of merging-cuts on the 
minimal energy of a cluster, E m i n , and on the separation distance between 
a pair of clusters, d pair . The merging-cut d pair is expressed in units of the 
Moliere radius, Rj^. The relative errors of the numbers for miscounted show- 
ers, iVi-^ and N2-+1, are computed for a sample of 3 • 10 4 Bhabha events. 
The number N to t takes into account both of the miscounting errors, and is 
computed for an integrated luminosity of 500 fb _1 . 
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shower (cluster) and multiplying by a calibration constant. The constant 
transforms between the values of the detector signal and the particle energy, 
and is determined by a calibration curve, such as the one shown in Fig. 5.6 
(Sect. 5.2.2). Figure 4.7 shows the normalized difference between the energy 
of reconstructed clusters and their respective generated showers, as a func- 
tion of the energy of the generated shower. The fluctuations are of 0(2 GeV) 
or lower. The reconstructed clusters, which are taken into account here, 
belong to the effective Bhabha cross-section for which E\ ow > 20 GeV and 
dpair > Rm- Increasing the merging-cut on separation distance reduces the 
fluctuations significantly, due to reduced cluster-mixing (see Sect. A. 3. 2), but 
also reduces the available statistics. 



e 0.2 r 

C3 

S! 0.15 r 

I 0.1- 
LU 

' E 0.05- 

-0.05 r 4 

-0.1 r 
-0.15 r 

■0 2 1 

U T) 50 100 150 200 250 

Ecen [GeV] 



Figure 4.7: The normalized difference between the energy of generated show- 
ers, E gen , and their respective reconstructed clusters, E rec , as a function of 
the energy of the generated shower. 



Figure 4.8 shows the normalized difference between the position of re- 
constructed clusters and their respective generated showers. The position is 
parameterized by the polar angle, 9, and the azimuthal angle, <fi. The differ- 
ence is presented as a function of the energy of the generated shower. The 
angles 9 and <fi are reconstructed according to Eqs. (A.l) and (A. 2). Since 
the fluctuations in all cases are of O(10 -5 ) or lower, it is concluded that 
the position reconstruction is performed well. This makes sense in light of 
Eq. (A.l), since only the core of high energy cells, which are in close proximity 
to the cluster center, contribute to the position reconstruction. Low-energy 
cells which are miss-assigned between clusters, therefore, do not degrade the 
reconstruction. 
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Figure 4.8: The normalized difference between the position of generated 
showers and their respective reconstructed clusters. The position is param- 
eterized by the reconstructed and the generated polar angles, 6 Rec and Gen , 
(a) and by the reconstructed and the generated azimuthal angles, 4>R ec and 
</>Gen, (b) and are presented as a function of the energy of the generated 
shower, E Gen . 
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Measurable Distributions 

The distribution of the energy of reconstructed clusters and their respective 
generated showers for high v and low-energy clusters (showers) is shown in 
Figs. 4.9a and 4.9b, respectively. 
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Figure 4.9: Distributions of the energy of reconstructed clusters (REC) and 
their respective generated showers (GEN), as denoted in the figures. The 
sample is divided into high (a) and low-energy (b) clusters (showers). 



Figure 4.10 shows the distributions of the polar angle, 6, of reconstructed 
clusters and their respective generated showers. The sample is divided into 
high and low-energy clusters (showers). In Fig. ?? is presented the distribu- 
tion of the difference in polar angle, A9 high j ow = 6 high — 9i ow , between the 
high and low reconstructed clusters and their respective generated showers. 

In light of the relations shown in Fig. 4.8 one might naively expect that 
the match between the distributions of cluster and shower positions would be 
better. The small noticeable discrepancies originate from miscounted events. 
On by-case basis the difference in reconstruction of the polar and az- 

imuthal angles is usually below the resolving power of LumiCal. However, 

v It may be noticed that the distribution of high-energy clusters has values smaller 
than 125 GeV (less than half the maximal possible energy). This is due to the fact that 
occasionally part of the energy is not included in the reconstruction. This may occur 
when some of the particles do not enter LumiCal, or when the position of a shower is 
reconstructed outside the fiducial volume, and so all cell-information is discarded. 
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Figure 4.10: Distributions of the polar angle, 9, of reconstructed clusters 
(REC) and their respective generated showers (GEN), as denoted in the 
figures. The sample is divided into high (a) and low-energy (b) clusters 
(showers). 




Figure 4.11: Distributions of the difference in polar angle, A8hi g h,iow, between 
the high and low reconstructed clusters (REC) and their respective generated 
showers (GEN), as denoted in the figure. 
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single showers which are reconstructed as two clusters, and shower pairs that 
are reconstructed as single clusters, must also be taken into account. The 
distributions in Fig. 4.6 indicate that the sources of the discrepancies are 
showers with small angular separation. This is indeed the case, as can be 
deduced from Fig. 4.12, where the instances of failure of the algorithm are 
shown as a function of the separation distance between pairs of showers, 
dp a i r , and of the energy of the low-energy shower, Ei ow . The merging-cuts 
Ei ow > 20 GeV and d pa i r > Rm have been used for selection of reconstructed 
clusters. The algorithm tends to produce mistakes when merging showers 
for which E\ ow and d pa i r are close to the merging-cut values, which is due to 
errors in either the position or the energy reconstruction. Thus, it is possi- 
ble to improve the results of the comparison between generated showers to 
reconstructed clusters, by making a selection-cut on events with low-cluster 
energy and cluster-pair distance, which are close to the merging-cut values. 
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Figure 4.12: Instances of failure of the clustering algorithm in distinguishing 
between a pair of generated showers, as a function of the separation distance 
between the pair, d pair , and of the energy of the low-energy shower, Ei ow . 
The distance d pa i r is expressed in units of mm and of Moliere radius, Rj^. 
Two cases are possible, a pair of generated showers may be merged into one 
reconstructed cluster (2 — > 1), or one shower may be separated into two 
clusters (1 — > 2). The event sample considered complies with the merging- 
cuts E low > 20 GeV and d pair > R M . 



4.4 Conclusions on Clustering 



47 



4.4 Conclusions on Clustering 



It has been shown that clustering of EM showers in LumiCal is possible. 
In order to achieve results of high acceptance and purity, a merging-cut on 
minimal energy for each cluster, and on the separation distance between any 
pair of clusters, needs to be made. The merging leads to a measurement of an 
effective Bhabha cross-section. The number of effective photons may then be 
counted with an uncertainty that corresponds to the required precision for the 
measurement of the luminosity spectrum. The distributions of the position 
and energy of the effective leptons and photons may also be measured and 
compared to the expected results. A merging-cut should be used in this case 
for summation of clusters, as is done for the counting of effective photons. 
Imposing an additional selection-cut on events can improve the results, by 
restricting even further the separation distance between the pair of clusters 
and the energy of the low-energy cluster, and thus effectively discarding most 
of the miscounted clusters. 



Chapter 5 

The Performance of LumiCal 



The performance of LumiCal may be evaluated using several parameters; 
the precision with which luminosity is measured, the energy resolution, the 
ability to separate multiple showers, viability of the electronics readout, and 
finally, the integration of LumiCal in the detector. In the following, each 
of these criteria will be discussed. The baseline geometrical parameters of 
LumiCal are presented in Table 5.1. It will be shown that for the present de- 
tector concept, these parameters fulfill the requirement of best performance 
of LumiCal. Finally, the influence of making changes to the different param- 
eters on the calorimeter performance will be summarized. This is necessary 
in order to facilitate setting an optimization procedure for future changes in 
the design of LumiCal. 



5.1 Intrinsic Parameters 
5.1.1 Energy Resolution 

LumiCal is designed in such a way that incident high energy electrons and 
photons deposit practically all of their energy in the detector. Energy degra- 
dation is achieved by the creation of EM showers, due to the passage of 
particles in the layers of tungsten (see Sect. 3.1). 

Prevention of leakage through the edges of LumiCal is possible by defining 
fiducial cuts on the minimal and on the maximal reconstructed polar angle, 
8, of the particle showering in LumiCal. Stable energy resolution is the 
hallmark of well-contained showers. The relative energy resolution, cte/E, is 
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Parameter 


Value 


Distance from the IP 


2270 mm 


Number of Radial divisions 


64 (0.75 mrad pitch) 


Number of Azimuthal divisions 


48 (131 rad pitch) 


Number of Layers 


30 


Tungsten thickness 


3.5 mm 


Silicon thickness 


0.3 mm 


Support thickness 


0.6 mm 


Layer gap 


0.1 mm 


Inner radius 


80 mm 


Outer radius 


190 mm 



Table 5.1: Baseline properties of LumiCal. 



usually parameterized as 



E 



E r 



beam 

(GeV) 



(5.1) 



where E and oe are the most probable value, and the root-mean-square of 
the signal distribution for a beam of electrons of energy Eb eam . Very often the 
parameter E res is quoted as resolution, a convention which will be followed 
in the analysis presented here. 

Figure 5.1a shows the energy resolution as a function of Q m in- The 
maximal angle is kept constant. The best energy resolution is achieved 
for Q m in = 41 mrad. A similar evaluation was done for a constant Q m %n 
and a changing 6 max , resulting in an optimal cut at 6 max = 69 mrad, as 
shown in Fig. 5.1b. The fiducial volume of LumiCal is thus defined to be 
within the polar angular range: 41 < 6 < 69 mrad. For this fiducial volume 
E res = (20.50 ± 0.05) • 10~ 2 y/(GeV). 
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Figure 5.1: Energy resolution for 250 GeV electrons as a function of the 
minimal polar angle, 8 m i n , (a) and as a function of the maximal polar angle, 



e, 



(b). 



5.1.2 Polar Angle Resolution and Bias 

The polar angle is reconstructed by averaging over the individual cells hit in 
the detector, using the cell centers and a weight function, Wj, such that 

< e >= ^» ' (5.2) 

Weights are determined by the so-called logarithmic weighting [7], for 
which 



W, = max{ , C + In— — }, (5.3) 

E-tot 

where Ei is the individual cell energy, E tot is the total energy in all cells, and 
C is a constant. In this way, an effective cutoff is introduced on individual 
hits, and only cells which contain a high percentage of the event energy 
contribute to the reconstruction. This cut, which depends on the size of the 
different cells, and on the total absorbed energy, is determined by C. There 
is an optimal value for C, for which the polar resolution, a$, is minimal. This 
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is shown in Fig. 5.2a using 250 GeV electron showers. The corresponding 
polar bias, A6, is presented in Fig. 5.2b. Accordingly, the polar resolution 
and bias of LumiCal are a e = (2.18 ± 0.01) • 10~ 2 mrad and A9 = (3.2 ± 
0.1) • 10~ 3 mrad. 
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Figure 5.2: The polar resolution, erg, (a) and the polar bias, AO, (b) as a 
function of the logarithmic weighing constant, C, using 250 GeV electron 
showers. 



5.2 Readout Scheme 

Upon deciding on the granularity of LumiCal, it is necessary to define the 
dynamical range of the electronics required to process the signal from the 
detector. Once the dynamical range is set, the digitization scheme depends on 
the ADC precision. The energy resolution and the polar-angle reconstruction 
depend on the digitization scheme. For the present study, it is assumed that 
the dynamical range of the electronics has to be such, that it enables to 
measure signals from minimum ionizing particles (MIP) up to the highest- 
energy EM showers, which are allowed by kinematics. 

In order to determine the lower bound on the signal in LumiCal, the 
passage of muons through the detector was simulated. Muons do not shower, 
and are, therefore, MIPs. In the present conceptual approach, muons will 
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be used to inter-calibrate the cells of the detector, and may also be used to 
check in-situ the alignment of the detector. The detection of muons in the 
forward region also has significance for many searches for physics beyond the 
Standard Model, such as implied by certain supersymmetry models, or by 
theories with universal extra dimensions [35]. 

In order to measure the signals of both MIPs and high energy electrons 
in LumiCal, the detector would have to operate in two different modes. In 
the calibration (high gain) mode the electronics will be sensitive to MIP 
signals. In the physics (low gain) mode the signals of high energy showers 
will be processed. The signature of a Bhabha event is an e + e~ pair, where 
the leptons are back to back and carry almost all of the initial energy. For 
the case of a nominal center of mass energy of 500 GeV, the maximal energy 
to be absorbed in LumiCal is, therefore, 250 GeV, and so 250 GeV electrons 
were used in order to find the upper bound on the detector signal. The low 
limit on the signal will have to be of the order of a single MIP, and will 
be precisely determined according to the restrictions imposed by the energy 
resolution. 

The output of Mokka is in terms of energy lost in the active material, 
silicon in the case of LumiCal. In order to translate the energy signal into 
units of charge, the following formula was used: 

1.6 ■ 10~ 4 

S [fC] = 3m § £ [eV] (5.4) 

where §£ denotes the signal in units of eV, and Sg the signal in units of fC. 
The value 3.67 eV is the energy to create an electron-hole pair in silicon. The 
number 1.6 • 10~ 4 fC is the charge of an electron. 

The detector model described in Table 5.1 was simulated. 
5.2.1 Dynamical Range of the Signal 

The distribution of the energy deposited in a detector cell by 250 GeV muons 
is presented in Fig. 5.3. According to this, the most probable value of induced 
charge for a muon traversing 300 /jm of silicon is 89 keV, which is equivalent 
to 3.9 fC. 

The distribution of collected charge per cell for 250 GeV electron showers 
is presented in Fig. 5.4a for a LumiCal with 96 or 64 radial divisions, which 
correspond to angular cell sizes of 0.5 and 0.8 mrad respectively. For the 
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Figure 5.3: Distribution of the energy deposited in a detector cell, E ce u, by 
250 GeV muons. A corresponding scale in units of charge is also shown. 



baseline case of 64 radial divisions, the value of the collected charge extends 
up to 6 pC, which is equivalent to ~ 1540 MIPs. The distribution of the 
maximal charge collected in a single cell per shower for 250 GeV electrons is 
shown in Fig. 5.4b for the two granularity options. As expected, for the case 
of smaller cell sizes the signal per cell is lower. 



5.2.2 Digitization 

Once a low or high bound on the dynamical range for each mode of operation 
is set, it is necessary to digitize the signal. For each mode of operation 
separately 



_ Hmax /- r\ 



I ma. r 



where ctadc is the ADC channel resolution (bin size), q m i n and q Tl 
respectively, the low and high charge bounds and 2' Bdi ^ is the number of 
channels for a given number of available ADC bits, "Bdigi- Each cell of de- 
posited charge, qd ep , is read-out as having a charge, qdigi, (rounding error) 
where 



qdigi = (Q(qde P , oadc) + 0.5) • oadc , 



(5.6) 
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Figure 5.4: (a) Normalized distribution of the charge deposited in a detector 
cell by 250 GeV electron showers for a LumiCal with 96 or 64 radial divi- 
sions, as denoted in the figure, (b) Normalized distribution of the maximal 
charge collected in a single cell per shower for 250 GeV electron showers for 
a LumiCal with 96 or 64 radial divisions, as denoted in the figure. In both 
figures a corresponding scale in units of MIPs is also shown. 
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and the quotient of the deposited charge with the ADC resolution, Q(qd ep , &xdc), 
is defined such that 



a = Q(a, 0) ■ (3 + 7 , 

0<7<|/3|- { } 

Table 5.2 shows the restrictions on the dynamical range of the two modes 
of operation. Since in the calibration mode the spectrum of MIPs will be 
measured, the resolution must be better than one MIP. For the choice of 
a low charge bound, q m in = 0.2 MIPs, the high bound will be determined 
by the digitization constant. The digitization constant will also determine 
the low bound for the physics mode, once the upper bound is set to q max = 
1540 MIPs. In Table 5.3 are presented the values of q max for the calibration 
mode and of q m i n of the physics mode for several choices of the digitization 
constant. 

Figure 5.5 shows the distribution of the total event energy for the case of 
a non-digitized, and an 8 or 10 bit digitized detector signal. The mean value 
of the distribution for the 8 bit digitized case is higher by 0(1%) compared 
to the non-digitized case. This is due to the fact that by far the greatest 
percentage of contributing cells have small energy, as can be observed in 
Fig. 5.4. The charge of more hits will, therefore, be overestimated, rather 
than underestimated. To illustrate this point for the 8 bit case, for the 
first ADC bin, 44% of the hits in a 250 GeV electron shower belong to the 
bottom half of the bin, while 21% belong to its top half. Since all of these 
hits are read-out with a digitized charge of exactly half the bin width, more 
contributions are overestimated, rather than underestimated, with weights 
according to the distribution of Fig. 5.3. One, therefore, finds that the total 
digitized charge is higher by 2.2 -r^— than the deposited charge, which 
amounts to a 26% increase in the integrated measured signal in this ADC 
channel. Since this effect depends on <tadc the difference in the mean value of 
the distributions of Fig. 5.5 with respect to the non-digitized case decreases 
for larger values of 'Bdigi- 

In practice, the shift in the mean bears no consequence other than the 
need to adjust the signal to energy calibration. The dependence of the detec- 
tor signal on the energy of the particle which initiated the shower is shown 
in Fig. 5.6 for several digitization schemes. There is no significant change as 
a result of the digitization, for the values of "Bdigi which were used. 

The important quantity that has to be controlled is the energy resolution, 
which must be the same for the digitized and the non-digitized cases. Fig. 5.7 
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Calibration Mode Physics Mode 

q min 0.8 fC (0.2 MIPs) a ADC 

q max 0.8 fC x2 1 W 6 pC (1540 MIPs) 



Table 5.2: Low and high bounds, q m i n and q m ax, of the dynamical ranges of 
LumiCal for operation in the calibration (high gain) and in the physics (low 
gain) modes. 



T> digi [bits] 


qmax of 

Calibration Mode 


qmin of 

Physics Mode 


6 


49.9 fC (13 MIPs) 


93.7 fC (24 MIPs) 


8 


199.7 fC (52 MIPs) 


23.4 fC (6 MIPs) 


10 


798.7 fC (205 MIPs) 


5.9 fC (1.5 MIPs) 


12 


3.2 pC (819 MIPs) 


1.5 fC (0.4 MIPs) 


14 


12.8 pC (3277 MIPs) 


0.4 fC (0.1 MIPs) 



Table 5.3: Dependence of the high bound, q m axi of the calibration mode, and 
of the low bound, q m i n , of the physics mode on the digitization constant, 
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Figure 5.5: Normalized distribution of the total event energy, E tot , of 250 GeV 
electron showers for a non-digitized, and either an 8 (a) or a 10 bit (b) 
digitized detector signal, as denoted in the figures. 
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Figure 5.6: Dependence of the detector signal, Etot, on the energy of the 
electron which initiated the shower, Ec e n- The detector signal is either un- 
digitized, or digitized with an 8 or a 10 bit digitization constant, as denoted 
in the figure. 
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shows the dependence of the energy resolution, E res , the polar resolution, a$, 
and the polar bias, A8 on the digitization constant. The values shown for 
^digi = 14 are equivalent to a non-digitized readout. 
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Figure 5.7: Dependence of the energy resolution, E res , (a) the polar resolu- 
tion, erg, (b) and the polar bias, A9, (c) on the digitization constant, "Bdigi- 



For CB d j 3 j > 8 it is apparent that the energy resolution, the polar res- 
olution, and the polar bias are all stable. Below this limit there is severe 
degradation of E res and slight improvement of oq and A6. The reason for 
this is that E res depends on the accuracy with which each and every detec- 



5.3 Geometry Optimization 



59 



tor cell is read-out. This means that fluctuations in q^gi of cells with small 
energy become critical when oadc is large in comparison to the cell signals. 
This effect, naturally, depends on the sizes of the LumiCal cells, since smaller 
cells have smaller signals, as suggested by Fig. 5.4. On the other hand, the 
polar angle reconstruction only takes into account contributions from cells 
with relatively large energy (Eqs. (5.2) and (5.3)), for which <tadc is small in 
comparison. Fluctuations, which hinder the polar reconstruction, decrease 
for low energy hits. With the negative influence on the energy resolution be- 
ing the driving factor, it is concluded that the minimization of E res requires 
the digitization constant to be higher than 8 bit. 

5.3 Geometry Optimization 

The chosen geometrical parameters for the baseline model are given in Ta- 
ble 5.1. The following is a systematic study, in which it will be shown that 
the values given in the table optimize the performance of LumiCal. Differ- 
ent single parameters will be varied, keeping the others constant, and the 
consequences of each change will be discussed. 

5.3.1 The Number of Radial Divisions 

For different radial granularity one needs to re-optimize the logarithmic 
weighing constant, C, of Eq. (5.3), as the size of cells changes for each case. 
The polar resolution and bias are plotted in Fig. 5.8 as a function of the 
angular cell size, £$. Their values are presented in Table 5.4, along with 
the corresponding relative error in the luminosity measurement (according 
to Eq. (2.9)). 

Both oq and A6 become smaller as the angular cell size decreases. The 
relative error in luminosity follows the same trend. This is due to the fact 
that the bounds on the fiducial volume do not strongly depend on the number 
of radial divisions. Consequently the minimal polar angle, 6 m i n , is the same 
(41 mrad) for all the entries of Table 5.4. 

Many problems arise when one increases the number of channels beyond 
a certain density. One has to resolve such problems as cross-talk between 
channels, power consumption issues, the need for cooling, etc. It is, therefore, 
advisable to keep the number of cells as low as possible. The other important 
parameter, the energy resolution, does not depend on the number of channels, 
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Figure 5.8: The polar resolution, erg, (a) and the polar bias, AO, (b) for the 
optimal logarithmic weighing constant, as a function of the angular cell size, 
in. Electron showers of 250 GeV were used. 



Radial Divisions 


ig [mrad] 


oq [mrad] 


AO [mrad] 


2A9 


96 


0.5 


1.8 ■ 10~ 2 


1.3 ■ 10~ 3 


0.6 • 10~ 4 


64 


0.8 


2.2 • 10" 2 


3.2 ■ 10" 3 


1.5 • 10" 4 


48 


1 


2.7- 10~ 2 


6.9 ■ 10~ 3 


3.1 • 10" 4 


32 


1.5 


3.5- 10~ 2 


13.7- 10- 3 


6.2 • 10~ 4 


24 


2 


4.4- 10~ 2 


24 ■ 10~ 3 


10.9 • 10" 4 


16 


2.5 


6.4- 10~ 2 


44.4 • 10~ 3 


20.2 • 10" 4 



Table 5.4: The polar resolution, erg, and bias, AO, for LumiCal with different 
numbers of radial divisions, corresponding to different angular cell sizes, ig. 
The corresponding values of the relative error in the measurement of the 
luminosity are also shown. 
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since the energy is integrated over all cells . The chosen baseline number of 
64 radial divisions is, therefore, a compromise between trying to minimize 
the relative luminosity error, and limiting the number of channels. 



5.3.2 The Structure of Layers 

Each layer of LumiCal consists of 3.5 mm of tungsten, which is equivalent 
to one radiation length. A distribution of the energy deposited in a layer by 
250 GeV electrons for a LumiCal of 90 layers is presented in Fig. 5.9a. Only 
0.4% of the event energy is deposited beyond 30 layers. The distribution 
of the total event energy for 250 GeV electrons is plotted in Fig. 5.9b for a 
LumiCal with 90 layers and for a LumiCal with 30 layers. A small difference 
is apparent in the mean value of the distributions, but this bears no conse- 
quence, as was explained in Sect. 5.2.2. The energy resolution is the same 
for both cases. Since there is no degradation of the energy resolution it is 
concluded that 30 layers are sufficient for shower containment. 




Layer E, ot [GeV] 

(a) (b) 

Figure 5.9: (a) Normalized distribution of the energy deposited in LumiCal as 
a function of the layer number for a detector with 90 layers, (b) Comparison 
of the distribution of the total deposited energy, E tot , for a 90 layer LumiCal 
with that of a 30 layer LumiCal, as denoted in the figure. 



This statement is true provided that the digitization constant is not too small, as 
discussed in Sect. 5.2.2. 
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The Moliere radius of LumiCal, Rm, is plotted in Fig. 5.10a as a function 
of the gap between tungsten layers. Since a smaller Ry& improves both the 
shower containment and the ability to separate multiple showers, the air gap 
should be made as small as possible. Fig. 5.10b shows the dependence of 
on the tungsten thickness. It is apparent that there is no significant change 
in i?M over the considered range. 
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Figure 5.10: The Moliere radius of LumiCal, Rm, as a function of the gap 
between tungsten layers (a), and as a function of the thickness of each layer 
(b). 



Changing the thickness of tungsten layers increases the sampling rate. In 
order to ensure shower containment, the total number of layers must remain 
30 radiation lengths. Consequently, for a thinner tungsten layer length, more 
layers are needed, as shown in Table 5.5. 



Tungsten Thickness [mm] 


2.0 


2.5 


3.0 


3.5 


4 


Number of Required of Layers 


53 


42 


35 


30 


26 



Table 5.5: The required number of LumiCal layers as a function of the thick- 
ness of each tungsten layer. 

Figure 5.11 shows the normalized distribution of the energy deposited per 
layer as a function of layer thicknesses, and the energy resolution, E res , for 
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each configuration. Figure 5.12 shows the corresponding polar resolution, erg, 
and bias, AO. 
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(a) (b) 

Figure 5.11: (a) Normalized distribution of the deposited energy for 250 GeV 
electron showers as a function of the layer number for several layer thick- 
nesses, as denoted in the figure, (b) The energy resolution, E res , as a function 
of the thickness of tungsten layers. 

Due to the fact that more layers encompass the shower-peak area for 
thinner tungsten layers, the energy resolution is improved. The polar recon- 
struction and the Moliere radius are not affected. The trade-off for choosing a 
given thickness of tungsten, is then between an improvement in E res and the 
need to add more layers. Since increasing the number of layers also involves 
a raise in the cost of LumiCal, a clear lower bound on E res needs to be de- 
fined, so as to justify the additional expense. Currently, E res « 0.21 ^/(GeV) 
seems sufficient. 

5.3.3 Inner and Outer Radii 

The distance from the IP (2.27 m), the radial cell size (0.8 mrad) and the 
number of radial divisions (64) dictate that the total radial size of LumiCal 
be 110 mm. Setting the inner and outer radii R m in and Rmax, within this 
limit has several implications. 
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Figure 5.12: The polar resolution, o~g, (a) and bias, AO, (b) as a function of 
the thickness of a LumiCal layer, using 250 GeV electron showers. 



In the detector integrated dipole (DID) [36] field configuration, the mag- 
netic field is directed along the incoming beam lines with a kink at the 
transverse plane containing the IP. Conversely, the magnetic field may also 
be directed along outgoing beam lines with a kink at the IP plane, a config- 
uration referred to as anti-DID. Figure 5.13 shows a projection of the energy 
of beamstrahlung pairs on the face of LumiCal for the anti-DID and the DID 
magnetic field configurations of the accelerator. The two inner concentric 
black circles represent possible inner radii of 60 and 80 mm, while the outer 
circle is set at 190 mm. The beamstrahlung spectrum was generated using 
GUINEA-PIG [37]. 

For a DID field the beamstrahlung pair distribution grazes LumiCal, while 
for the anti-DID case it does not, though the distribution comes close to the 
LumiCal inner edge. The affect of exposure to the extremely high energy 
dose will cause massive damage to the silicon sensors in a matter of months. 
The anti-DID field is, thus, the better choice. The difference between the 
two distributions of Fig. 5.13 suggests that small fluctuations in the magnetic 
field from the nominal configuration will cause the pair distribution to become 
wider. It is, therefore, concluded that it would be preferable to add a safety 
margin to the minimal (60 mm) choice of R m in- 

It has also been shown [38] that for R m i n < 70 mm there is a significant 
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Figure 5.13: Projection of beamstrahlung pair energies on the face of Lumi- 
Cal for an anti-DID (a) and a DID (b) magnetic field setup, for the nominal 
accelerator operational parameters. The concentric circles represent possible 
inner LumiCal radii of 60 and 80 mm, and an outer radius of 190 mm. 



increase in the amount of backscattered particles from LumiCal to the inner 
detector (the TPC). This too constitutes a motivation for setting R m i n at a 
higher value. 

Since the Bhabha cross-section falls off quickly with the polar angle (see 
Eq. (2.4)), it is advantageous to set R m in as low as possible in order to 
increase the number of Bhabha events within the fiducial volume of LumiCal. 
Table 5.6 gives the integrated Bhabha cross-section, <tb, in the fiducial volume 
defined by several choices of R m in and R max - The number of Bhabha events 
and the relative statistical error is calculated according to Eqs. (2.1) and 
(2.10), respectively. An integrated luminosity of 500 fb~ was assumed. The 
relative error resulting from the polar reconstruction (Eq. (2.9)) is also shown 
in the table, where a polar bias A6 = 3.2 • 10~ 3 mrad (Table 5.4) and the 
appropriate minimal polar angles were used in each case. 

As expected from Eq. (2.4), for low values of R m in the number of Bhabha 
events increases, thus decreasing the statistical error. As the polar bias 
depends on the angular cell size, which was kept constant, and not on the 
radii, the error resulting from the polar reconstruction decreases slightly with 
the rise of Q m in- Both effects contribute to the overall uncertainty in the 
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Rmin ^ Rmax 


@min 


@max 




AN B 


2A6> 




[mm] 


[mrad] 


[mrad] 


[nb] 


N b 






60 ->■ 170 


33 


59 


2.58 


2.8 ■ 10~ 5 


1.9-10" 


-4 


70 ->■ 180 


37 


64 


1.98 


3.2 ■ 10~ 5 


1.7 - 10 


-4 


80 ->■ 190 


41 


69 


1.23 


4 • 10~ 5 


1.5 ■ 10 


-4 


90 ->■ 200 


50 


74 


0.86 


4.8 ■ 10~ 5 


1.3-10" 


-4 



Table 5.6: The fiducial volume, R m in — ► Rmax-, bound by the minimal and 
maximal polar angles, Q m in and 8 max , and the integrated Bhabha cross- 
section, Ob, for a center of mass energy y/s = 500 GeV. The two relative 
errors on the luminosity measurement, the statistical error and the one re- 
sulting from reconstruction of the polar angle, are also shown. The number 
of Bhabha events, Nb, is computed for an integrated luminosity of 500 fb" 1 , 
and the polar bias used is A6 = 3.2 • 10~ 3 mrad. 

luminosity measurement. 

It should be noted here that in practice the counting rates of Bhabha 
events will be lower than presented in Table 5.6. This is due to the fact 
that the efficiency for counting Bhabha events is not 100% due to selection 
cuts [25]. This in itself does not add to the luminosity error, as long as 
the efficiency is known to high precision, but it does increases the statistical 
error 11 . 

The contribution of the error due to the polar bias also needs further 
consideration. In practice it will be possible to determine the polar bias using 
a test beam, and correct for this effect. The final error will then depend on 
how well one can correct for the bias, so that the values given in the table 
are an upper bound on the error. 

In conclusion, even given the lower counting rate and no corrections of 
the polar bias, it is apparent that for the range of R m in — > Rmax given in 
Table 5.6, the relative error in luminosity is well within the design goal. Aim- 
ing to increase the available statistics as much as possible, while maintaining 
a safe distance from the beamstrahlung pairs, R m %n — 80 mm was finally 
chosen. 



11 For instance, for a pessimistic selection efficiency of roughly 50%, the error will increase 
by a factor of \/2. 
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5.3.4 Clustering 

The effectiveness of the clustering algorithm, which was described in chapter 4, 
in measuring the Bhabha cross-section, depends on the granularity of Lumi- 
Cal. The performance, which was presented in Sect. 4.3.2, was evaluated for 
a detector with inner and outer radii, R m in = 80 mm and R ma x = 350 mm 
respectively, 104 radial divisions, and 96 azimuthal divisions. This geometry 
corresponds to a radial cell size of 2.6 mm, and an azimuthal cell size of 
8.5 mm at the center of LumiCal. For the optimized geometry, described in 
Table 5.1, the radial cell size is 3.25 times smaller, and the azimuthal cell 
size is twice as large. 

In order to estimate the dependence of the performance of the clustering 
algorithm on these changes, the clustering of a sample of 10 4 Bhabha events 
was performed for different LumiCal geometries 111 . In all cases the inner and 
outer radii were kept at the optimized values of R m in = 80 mm and R ma x = 
190 mm, respectively, and the number of radial and azimuthal divisions was 
changed. The merging cuts used on the minimal energy of a cluster and 
on the separation distance between a pair of clusters (see Sect. 4.3.1) are 
Eiow > 20 GeV and d pa i r > Rm, respectively. 

The values of the acceptance and purity (Eq. (4.1)) are presented in 
Table 5.7. Also shown is the total relative error on the measurement of 
the effective Bhabha cross-section (Eqs. (4.3) and (4.4)) for an integrated 
luminosity of 500 fb _1 . 

Changes in the number of azimuthal divisions have a large effect on the 
final error of the cross- section measurement, compared to changes in the 
number of radial divisions. This difference is due to the fact that LumiCal 
is more finely granulated in the radial direction. One may fine-tune the 
parameters of the algorithm in order to adjust the values of the acceptance 
and of the purity. In general, an increase in A will be followed by a decrease 
of CP, as one is a measure of over-merging of clusters, and the other of under- 
merging. The contribution to the statistical error of the number of shower- 
pairs which are reconstructed as one cluster, far outweighs that of the number 
of single showers which are reconstructed as two clusters. Since it is advisable 
to choose parameters, such that the total error is minimal, the acceptance 
tends to be higher than the purity. 

111 The clustering algorithm described in Appendix A was utilized. Adjustments of several 
of the parameters were made, such as the number of near-neighbors and the weighting 
constant, C, (Eq. (A. 2)) in order to accommodate the changes in cell size. 
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Number of Divisions 


A [%] 


7 [%} 


AN tot 


Azimuthal 


Radial 


Jv tot 


96 


128 


99 


94 


2.9 ■ 10~ 5 


96 


32 


98 


92 


3.5- 10" 5 


48 


128 


94 


79 


6.6 ■ 10" 5 


48 


64 


93 


77 


7.5- 10" 5 


48 


32 


90 


84 


9.1 ■ 10" 5 


24 


64 


76 


22 


11.1 ■ 10~ 5 



Table 5.7: The values of the acceptance, A, and purity, CP, of the algorithm, 
for several division schemes. Also shown is the total relative error on the 
measurement of the effective Bhabha cross-section for an integrated lumi- 
nosity of 500 fb _1 . The merging-cuts on the minimal energy of a cluster and 
on the separation distance between a pair of clusters are E[ ow > 20 GeV and 
dpair > Rm, respectively. 

It is apparent that the increase in cell size diminishes from the effective- 
ness of the clustering. For the optimized geometry, though, the final relative 
error with which the Bhabha cross-section may be measured, is within the 
design goal. 

5.4 Conclusions on the Optimization Proce- 
dure 

In the following, the dependence of the various performance parameters on 
the geometry of LumiCal are summarized. 

Dynamical range of the signal - Reducing the maximal signal of a single 
cell may be accomplished by increasing the number of cells (Fig. 5.4). A 
low signal size brings about less rounding errors when digitizing the induced 
charge in a cell, and is also preferable from the readout electronics point of 
view. This change, however, carries the added complication of increasing the 
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number of channels, which in turn hinders the readout, as discussed above. 

Energy resolution - The energy resolution depends on the containment of 
the EM shower, on the precision with which each cell is read-out, and on the 
sampling rate of the shower. 

- Containment of showers is achieved by keeping the total number of 
layers 30 X thick, and imposing fiducial cuts on the polar angle of 
incident showers. 

- The accuracy of reading-out cell-energies is guarantied not to degrade 
the energy resolution, so long as a digitization scheme with high enough 
resolution (Fig. 5.7a) is implemented. 

- The sampling rate of the shower is determined by the thickness of 
each tungsten layer (Fig. 5.11). The best way to improve the energy 
resolution is to decrease the thickness of layers. One must increase the 
number of layers accordingly. 

Moliere radius - Keeping the Moliere radius small improves both the shower 
containment, by relaxing the fiducial cuts, and the ability to resolve multiple 
showers. This may be done by decreasing the gap between tungsten layers 
to the minimal possible value. 

Inner and outer radii of LumiCal - The inner radius of LumiCal deter- 
mines the minimal polar angle, which is accessible by the calorimeter. Since 
the Bhabha event-rate falls off rapidly with the polar angle (Eq. (2.4)), it is 
preferable to decrease R m %n as much as possible. The lower bound on R m i n 
must be set such that the beamstrahlung pairs do not enter LumiCal, and 
the backscattering from LumiCal into the inner detector is acceptable. The 
outer radius of LumiCal is less important in terms of a gain in statistics. 

Relative error on the luminosity measurement - Possible future changes 
in the position and size of LumiCal may affect its angular coverage and its 
angular cell size. This will influence the accuracy of the polar angle recon- 
struction, and accordingly the relative error on the luminosity measurement 
(Fig. 5.8 and Table 5.4). The number of radial divisions must therefore be 
adjusted, so that the angular cell size remains constant. Decreasing the num- 
ber of azimuthal cells does not affect the reconstruction of the polar angle, 
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yet a measurement of the Bhabah cross-section through clustering, requires 
that cell sizes not be large in comparison to the Moliere radius. The total 
fiducial volume of LumiCal must also be large compared to Rm, so as to 
achieve a high rate of Bhabha events, which pass the merging-cut on the 
separation distance between cluster pairs. 



Chapter 6 



Summary 



The International Linear Collider will provide physicists a new doorway to 
explore energy regimes beyond the reach of today's accelerators. A proposed 
electron-positron collider, the ILC will complement the Large Hadron Col- 
lider, a proton-proton collider at the European Center for Nuclear Research 
in Geneva, Switzerland, together unlocking some of the deepest mysteries 
in the universe. With LHC discoveries pointing the way, the ILC, a true 
precision machine, will provide the missing pieces of the puzzle. In order 
to achieve its goal, the luminosity of the ILC will have to be know with a 
precision of 10 -4 , which poses a significant challenge. Luminosity in the ILC 
is measured by counting Bhabha scattering events, a task which will require 
pattern recognition in the main detector, and in the luminosity calorimeter. 

It has been shown here that it is possible to resolve the distribution 
of radiative Bhabha photons on top of the electron distribution, and thus 
measure an effective Bhabha cross-section directly. Using this measurement, 
it will be possible to verify the influence of the beam-beam effects, and of 
the energy spread of the collider, on the Bhabha cross- sect ion. 

A study has also been presented, in which the design of LumiCal was 
optimized, with the goal of reducing the error in the luminosity measurement 
to the required threshold. While this study holds merit in and of itself, it also 
serves as a template for future optimizations. Such procedures are foreseen 
to be needed, in light of the expected changes to the detector-concept, as a 
result of ongoing R&D efforts. 
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Appendix A 



The Clustering Algorithm in 
Full Detail 



The clustering algorithm, with was developed for LumiCal, operates in three 
main phases; 

- selection of shower-peak layers, and two-dimensional clustering therein; 

- fixing of the number of global (three-dimensional) clusters, and collec- 
tion of all hits onto these; 

- testing of the global-clusters, by means of the evaluation of their energy 
density. 

The algorithm is described in detail in the following subsections. 

A.l Clustering in the Shower- Peak Layers 
A. 1.1 Near-Neighbor Clustering 

As the initiating particle traverses LumiCal, more and more secondary par- 
ticles are created and the shower spreads out (see chapter 3). The greatest 
density of energy per individual layer is concentrated around a local center-of- 
gravity. The method of near-neighbor clustering (NNC) exploits the gradient 
of energy around the local shower- center, by assuming that in first order, the 
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further a hit is relative to the shower center, the lower its energy. By com- 
paring the energy distribution around the center at growing distances, one 
may check whether the energy is increasing or decreasing. An increase in 
energy for growing distance from the shower-center would then imply that 
the hit should be associated with a different shower. 

It is helpful to consider at this point only cells which have a relatively 
high energy content, and so a cut on cell energy (given below in Eq. (A. 3)) 
is made before the NNC begins. Hits of low energy are, therefore, associated 
with clusters at a later stage. 

The radial cell length is constant and has a value of 2.6 mm (for the geom- 
etry described in Sect. 2.2), while the azimuthal cell length, which is radially 
dependent, is 5 mm at its lowest and 23 mm at its highest. The Moliere ra- 
dius of LumiCal is Rm = 14 mm, which is comparable with the cell lengths in 
the middle of the detector. This means that the cell size of LumiCal hinders 
the effectiveness of the procedure described above, and therefore the NNC 
method is used in a mode, in which only six near-neighbors are considered 
for each hit. One and two displacements are considered in the radial direc- 
tion, and one in the azimuthal direction. A schematic representation of the 
relative locations of the near-neighbors is shown in Fig. A.l. The lighter the 
shade of the cell, the further away is its center, relative to the center of the 
principle cell. The ratio of the radial to azimuthal cell size is representative 
of the middle region of LumiCal. 

For each shower-peak layer separately, the algorithm associates each cell 
which has an energy deposit with its highest-energy near-neighbor. The 
result of the NNC phase is a collection of clusters in each layer, centered 
around local maxima, as illustrated for a single layer in Fig. 4.1a (Sect. 4.1.1). 
In this example the algorithm produces six clusters, which are enumerated 
in the figure. The different clusters are also distinguished by different color 
groups, where darker shadings indicate a higher energy content of the cell in 
question. 

A. 1.2 Cluster-Merging 

The next phase of the algorithm is cluster-merging. The NNC method only 
connects cells which are relatively close, while showers tend to spread out 
over a large range of cells, as indicated in Fig. 3.5b (Sect. 3.2). The result of 
the algorithm after the cluster-merging phase (explained below) is illustrated 
for a single layer in Fig. 4.1b (Sect. 4.1.1). 
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Figure A.l: Schematic representation of the relative locations of the near- 
neighbors which are considered in the NNC algorithm. Darker shadings 
represent shorter distances between the cells in question and the principle 
cell, which is at the center of the drawing, painted black. 



The first step in the cluster-merging phase is the assignment of a center- 
position for each cluster. This is done by averaging over all the hits of each 
cluster, using the hit cell centers, (xi,yi), and a weight function, Wj. For the 
x-coordinate of the shower-center one has 



<I>= ^wT' (A1) 

and equivalently for the y-coordinate. Weights are determined by the so- 
called logarithmic weighting [7], for which 



Wi = max{ , C + In— -}, (A.2) 

E c i 

where is the individual cell energy, E c i is the sum of the energy of all 
the cluster cells, and C is a constant. The constant, C, acts as an effective 
cut on energy, so that only cells which contain a high percentage of the 
cluster energy contribute to the reconstruction of the position of the cluster. 
Setting different values for C influences the performance of the algorithm. It 
was found that the value of C, determined for a single shower of 250 GeV 
(Sect. 5.1.2) also gives the best results for the case of multiple showers with 
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the same total energy. It was also found that one should use the same 
weighting constant for all of the showers. 

Inverting Eq. (A. 2) allows to finally define the cut on cell energy, C m i n , 
which was used before the NNC began, 

Cmin{7 cut) = G ^ E a u ■ 7 cut > (A. 3) 

where E a u is the total energy deposited in all the layers of LumiCal, and 7 cu t 
is some fraction. The best performance of the algorithm was achieved for 

Once a center-position is determined for each cluster, the merging pro- 
cess begins. A vector of weights, W merge , is computed for each cluster, the 
elements of which are the weights, W^ erse , for the cluster, i, to merge with 
a cluster j (for i, j £ [1, n] , i ^ j). Weights are defined such that 



I (J otherwise. 

where Ej is the total energy of cluster j, and dij is the distance between 
the two clusters. The symbols a and (3 are parameters of the weighting 
process and ^l V cject is a rejection condition. Clusters are arranged according 
to their energy in an ascending order, and W merge is then computed for each 
cluster. A cluster is merged with the partner j for which W^ er9e £ W merge is 
maximal. In the case that {Vj £ n : W^ er9e = 0}, no merging is performed. 

The rejection condition is determined according to two criteria. The first 
is that the two candidates must not be too far apart, and the second is that 
the energy density of the merged cluster must be high, close to its cluster- 
center. Fig. 4.1b (Sect. 4.1.1) illustrates this point well. Had clusters one 
and two been mistakenly merged, the cluster-center would have ended up 
being in the space between the two, where the energy density is zero. The 
rejection condition may be written as 

Qreject(d m im ^sum) ^ 
if *C drain Or &cl{dsum) ^ E mer g ec i , 

where d m i n is the minimal separation distance between two clusters, E merge d 
is the total energy of the merged cluster, and t c i is the amount of energy 
contained within a distance d sum of the merged cluster-center. 



A.l Clustering in the Shower-Peak Layers 



77 



The best performance of the algorithm was achieved using: 



'mm 




(A.6) 



d 



sum 



Each cluster is defined as either small or large, according to the number 
of cells which the cluster holds, compared to the total number of hits in the 
layer. Small clusters are defined as having no more than 10% of the hits in the 
layer, and large clusters as having no less than 15%. Merging is attempted 
twice. In the first pass, small clusters attempt to merge with large clusters, 
and in the second, large clusters attempt to merge with the small. In this 
way the merging process becomes smoother, as the addition of large clusters 
to small tends to change the cluster-center of the small cluster much more 
than that of the large. It was also found that for the merging of small clusters 
with large, the distance between the pair is more important than the energy 
of the large cluster. For the case of the merging of large clusters to small 
the conclusion was that the two factors are equally important. The optima 
of values for the weighting constants, a and (3, for merging small clusters 
with large (1 st pass) and for merging large clusters with small (2 nd pass) are, 
therefore, 



A. 1.3 Clustering of the Low-Energy Cells 

Once the cluster-merging is over, the low-energy cells, which were discarded 
previously, are incorporated into existing clusters. This is done using weights 
computed according to the inverse of the distance between the hit in question 
and each of the clusters in the layer. Weights are, therefore, determined such 
that 




and 



(A.8) 



(A.7) 



-i 



(A.9) 



merge 
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which is equivalent to Eq. (A. 4) with parameters a = and (3 = — 1, and no 
rejection condition. 

A. 2 Global Clustering 

The basic actions taken in the global-clustering phase £1X6 cLS follows: 

1. fixing the number of global-clusters according to a frequency count in 
the shower-peak layers; 

2. re-clustering of hits from excess clusters where there are too many layer- 
clusters, and discarding of all clusters from layers where there are too 
few layer-clusters; 

3. creation of virtual-clusters (with a given center-position and size) in 
non-shower-peak layers, and in those where the layer-clusters were dis- 
carded; 

4. clustering of all remaining hits into the virtual-clusters, which are then 
identified with global-clusters. 

The most important stage of the clustering algorithm is the determination 
of the number of reconstructed showers. The aftermath of clustering in the 
shower-peak layers is several collections of two-dimensional hit aggregates, 
the number of which, n c i(£), varies from layer to layer. The final number of 
showers, N c i, is then determined as the most frequent value of n c i(£) from 
the collections in the shower-peak layers. 

The first step in the global-clustering procedure is calculating the center- 
position of the global-clusters. Since layers for which n c i{£) ^ N c i introduce 
ambiguities, they are temporarily not considered. In all accepted-layers, 
where n c i(£) = N c i, matching ensues between the layer-clusters and the 
global-clusters. This is done by comparing the two-dimensional distance 
between the centers of layer-cluster pairs in different layers, and choosing the 
minimally separated pairs. With all clusters matched, the position of the 
global-clusters is computed using Eqs. (A.l) and (A. 2), where E c i now refers 
to the total energy of the global-cluster. 

Turning back to the shower-peak layers for which n c i(£) ^ N c i, two sit- 
uations are possible. For layers where n c i(£) > N c i, the N c i best matches 
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are made between the layer clusters and the global clusters. The remaining 
(^d(^) — N c i) clusters are disbanded, and their hits are associated with one 
of the N c i remaining clusters. Once again, this is done using Eq. (A. 9). 

Regarding shower-peak layers where n c i(£) < N c i, all information gained 
by the algorithm is ignored. For these and for the non shower-peak layers 
clustering is performed in the following manner: firstly, shower-centers are 
fixed in all un-clustered layers. This is done by fitting a straight line through 
the centers of layer-clusters in the accepted layers, and extrapolating the 
position onto the rest. In this way, each of the non-accepted-layers now 
has N c i virtual- clusters, which correspond to the different global-clusters. 
Secondly, in each of the non-accepted layers separately, individual hits are 
associated with one of the virtual-clusters. 

Adding hits to the virtual-clusters is done in several steps. In the first 
stage, the virtual- cluster range is estimated to be the area around the virtual- 
center, which is spanned by the effective layer-radius of the layer, r e ff(£) 
(Eq. (3.2) in Sect. 3.2). All of the hits inside virtual-cluster ranges are added 
to their respective virtual-clusters. Cases where a hit is within the range of 
more than one virtual-cluster are resolved by a proximity test to the virtual- 
center positions. Once a core of hits has been added to each virtual-cluster, 
the center-position of each one is computed in the usual manner. All of 
the remaining un-clustered hits are finally associated with one of the now 
real-clusters by means of Eq. (A. 9). It is important to perform this process 
in these two phases, i.e., first to gather hits from the close proximity of 
the virtual-center, and only then to merge all of the remaining hits. The 
reason for this is that the real-cluster center, the position of which is mostly 
determined by its core hits, may be slightly different than the extrapolated 
virtual-center. This difference tends to bias the entire procedure in a given 
layer, if it is not accounted for. 

A schematic representation of the global-clustering phase is depicted in 
Fig. 4.2 (Sect. 4.1.2). 

A. 3 Corrections Based on the Energy Distri- 
bution 

By this point all of the hits in the calorimeter have been integrated into one 
of the global-clusters. Before moving on, it is beneficial to make sure that 
the clusters have the expected characteristics. 
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A. 3.1 Energy Density Test 

EM shower development in LumiCal has been described in Sect. 3.2. Ac- 
cordingly, one would expect that 90% of the energy of a cluster would be 
found within one Moliere radius, Rm, of its center. While statistically this is 
true, on a case-by-case basis fluctuations may happen, and thus this should 
not be taken as a hard rule. One can define, instead, the following general 
conditions: 

N c i N cl 

£ci{Rm) >J2 Ed - ' and ( A - 10 ) 

cl cl 

£d(i?M ■ Vci) > E d ■ ? tot , (A.ll) 

where Nd is the number of clusters in the layer, E c i{d) is the energy in a 
cluster, cl, within a distance d of its center, E c i is the total energy of cluster 
cl, and 7 tot and T c i ar e fractions. 

As an initial check, we evaluate Eq. (A. 10) with the constant y tot = 
0.9, meaning that the energy around each of the reconstructed clusters is 
integrated around Ryu of its center, and the sum of these is compared to 90% 
of the total energy found in the detector. If the condition is not upheld, 
an attempt is made to re-cluster and produce profile- clusters, based on the 
profile image of the energy distribution (see Fig. 3.4 in Sect. 3.2). The idea is 
that the core of the hits of a shower tends to stand-out against the background 
of low-energy deposits. By stripping the noise away, one is able to trace the 
propagation of the particle, which initiated the shower, through the different 
layers, and conduct a rough position reconstruction. The following steps are 
taken: 

1 . Hit energies from all layers are integrated along the longitudinal direc- 
tion, resulting in a two-dimensional energy distribution. 

2. The NNC method (Sect. 4.1.1) is utilized in order to build the profile- 
clusters using all of the available hits. 

3. The condition of Eq. (A.ll) is evaluated with the constants 7^ = 0.4 
and % ot = 0.8. 

4. If the condition of step (3) is not upheld, then a cut on hit-energy is 
made, removing some percentage of the lowest-energy hits. 
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5. Steps (2)-(4) are repeated until either the condition of step (3) is met, 
or the number of remaining hits after the progressive energy cuts drops 
below a certain threshold, in which case the construction of profile- 
clusters has failed. 

6. If the condition of step (3) is met, then all of the discounted low-energy 
hits are taken back into consideration, and Eq. (A. 11) is evaluated with 
^tot = 0.9. If this is upheld, the profile-clusters are accepted. 

If the building of profile-clusters is successful, the resulting shower-centers 
of the profile-clusters are used in order to construct global-clusters. This is 
done on a layer-to-layer basis, by means of the weights given in Eq. (A. 9). 
The new global-clusters are then compared to the original global-clusters by 
computing for each cluster set the energy density, pe, which is defined as 



If the profile-clusters have a higher energy density, then they are an improve- 
ment on the original clusters, and are kept instead of the originals. 

In addition to the procedure described above, each cluster is also exam- 
ined separately, with a relaxed condition (A. 11), using 7 c i = 1 and 7 tot = 0.5. 
Clusters that fail the test are disbanded, and their hits are assigned to other 
existing clusters using Eq. (A. 9). 



Another modification that can be made in the aftermath of the clustering 
procedure, is allocation of hits for mixed cluster pairs. A convenient way to 
define the intermixing of a pair of clusters, is to project the hits belonging to 
each cluster on the axis defined by the cluster-centers. The projection axis is 
schematically represented in Fig. A. 2. All of the hits which were represented 
in the original coordinate system, (x, y), are integrated along the y' direction 
of the projection coordinates, (x',y'). 

The distribution of the projected energy for a sample of cluster pairs is 
shown in Fig. A. 3. Several hundred events were simulated, and the projected 
energy distributions are all normalized according to the high-energy shower 
distribution, so as to keep the ratio between distributions constant. The 
energy of the initiating particles for the large and small-energy showers is 




(A.12) 



A. 3. 2 Unfolding of Mixed Clusters 
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230 and 20 GeV, respectively. Both showers were simulated according to a 
flat distribution along the detector face, while keeping the separation distance 
between the two at a constant d pa i r = 20 mm (~ 1.5 • Rm.) on the first layer. 
This way, the effects of different inclinations with regard to cell pitch, as well 
as the shower-shape distortions that are expected due to changing cell size, 
are smoothed out. In order to produce a clean distribution, the coordinates 
of the high-energy shower-center were always translated to (x',y') = (0,0), 
and those of the low-energy shower to a positive value along the projection 

clXIS, X . 

The distributions in Fig. A. 3 share a mixing-range, that extends from the 
peak of one shower to the next (0 < x' < d pair ). This, of course, is due to 
the fact that the un-projected distance between the pair is one and a half 
Moliere radius, beyond which the energy content of each shower is negligible. 
This fact poses an opportunity for resolving the mixing of the hits of the 
two showers. The area of negative x' holds a clean sample of hits, which 
have energy contributions from the high-energy shower alone. For x' > d pa i r , 
only the low-energy shower contributes significantly. It is, therefore, possible 
to estimate the amount of energy deposits owing to each shower inside the 
mixing-range, by distracting from the mixed distribution either one of the 
clean distributions, (x' < 0) or (x' > d pa i r ). 

For this to work, one has to assume that the shape of the original (un- 
projected) energy distributions is smooth. This is not always the case, 
though, as fluctuations in the development of showers tend to distort their 
shape. Other sources of distortions are the inclination of the projected axis, 
x', with relation to the cells of LumiCal, and the radial dependence of cell 
sizes, as mentioned above. 

The fact that the reference for corrections is the projected energy dis- 
tribution also poses problems, as one would like to know the original two- 
dimensional distribution of energy for each cluster. The proper assignment 
of individual cell energies, as opposed to a global correction on the inte- 
grated cluster energy, is important. The reason for this is that the amount 
of energy in each individual cell determines its weight in the position re- 
construction procedure (Eq. (A.2)). For the low-energy cluster especially, 
wrongly assigned cells of high energy, which comprise a significant amount of 
the shower energy, may bias the position of the reconstructed cluster-center. 

While the ambiguities mentioned here prevent from drawing global cor- 
rection factors directly from the distributions of Fig. A. 3, on a case-by-case 
basis, some degree of unfolding is possible. The way this is done is by project- 
ing the hits of every cluster pair in the (x', y') system, as described above, and 
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Figure A. 2: Schematic representation of the definition of the projection 
coordinate system, (x',y'), between a pair of intermixed clusters. 




x' [mm] 

Figure A. 3: Projected energy distribution for a sample of high and low- 
energy showers, initiated by electrons of 230 and 20 GeV, respectively. The 
distributions are normalized according to the high-energy shower distribu- 
tion. The two particles are evenly distributed on the face of LumiCal, and 
are generated with a constant separation distance of 20 mm on the first 
layer. The coordinates of the high-energy shower-center are translated to 
(V, y') = (0, 0), and those of the low-energy shower to a positive value along 
x', according to the separation between the two showers. 
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then determining weights according to the ratio of the energy distributions 
at x' > and x' < 0. These weights are defined as 



VWfc) = EproA f Jp^ , (A.13) 



where 



W&TH Y, V W oy (A-14) 

The weights, W mi;E (6j), are calculated for bins, &j, along the projection axis 
for (0 < < d pair ), where d pair is the distance between the centers of the 
two clusters. Energy is integrated from all hit contributions, £ pro j, in bin 
pairs of equal distance from the origin {x' = 0) along the positive (6f ) and 
negative (bf) directions. Bin sizes were chosen to be 5% of Rjy[. 

Returning to the un-projected energy distributions, each hit that would 
have satisfied (0 < x' < d pair ) had a projection been made, is split. A fraction 
of the hit energy, determined by W m j X (fej), is associated with each of the two 
clusters. The low-energy cluster looses W mix (6j)% of the energy of each hit, 
and the high-energy cluster gains the same amount. Splitting of hits is only 
performed for cases where *W m i x (bi) > 0, i.e., where more energy is projected 
in the mixing-range of the two clusters, than in the respective clean range of 
the high-energy cluster. 

The option of comparing the distribution of energy in the mixing-range 
to that of the clean range of the low-energy cluster (x 1 > d pa i r ) was also 
considered. This did not produce stable results, though, due to the fact 
that low-energy clusters tend to fluctuate more in shape, and are also more 
susceptible to cell geometry changes. 

Figure A. 4 shows the projected energy distributions of the generated 
(simulation-truth) and reconstructed (results of the clustering) energy dis- 
tributions before and after the unfolding of hits. The sample is the same as 
the one which was shown in Fig. A. 3, and is normalized in a similar manner. 
It is apparent that the splitting of hit energies improves the distribution of 
energy between the two clusters. It may also be noticed from the figure that 
the energy of low-energy clusters tends to be overestimated when mixing is in 
play, which explains why no splitting of hits is performed for W m j X (&i) < 0. 

The improvement of the energy reconstruction, due to the unfolding pro- 
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Figure A. 4: (a) Distributions of the energy in the projection coordinate 
system for the high-energy and low-energy showers/clusters. Both the 
simulation-truth shower (GEN) and the reconstructed cluster (REC) dis- 
tributions are shown. No unfolding of the distributions has been performed, 
(b) Distributions of the projected energy of the original (REC - Orig), and 
of the unfolded (REC - Fixed) low-energy reconstructed clusters. The dis- 
tributions of the simulation-truth (GEN) are also shown for the low and 
high-energy showers. In both (a) and (b) the distributions are normalized 
according to the high-energy shower distribution. 
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cedure, can be quantified by studying its affect on a physics sample of Bhabha 
events (see Sect. 4.2). Figure A. 5 shows the normalized difference between 
the energy of reconstructed clusters and their respective generated showers, 
as a function of the energy of the generated shower, before and after the un- 
folding procedure. The reconstructed clusters, which are taken into account 
here, belong to the effective Bhabha cross-section for which Ei ow > 20 GeV 
and dp a i r > Rm- 
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Figure A. 5: The normalized difference between the energy of generated show- 
ers (GEN) and their respective reconstructed clusters (REC), as a function 
of the energy of the generated shower, before (a) and after (b) the unfolding 
procedure. 
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